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ABSTRACT — In this paper (in Ralian) we critically and detaily examine various def- 
oo : initions existing in the literature for the tunnelling times: namely, the phase-time; the centroid- 
based times; the Biittiker and Landauer times; the Larmor times; the complex (path-integral 
and Bohm) times; the dwell time, and the generalized (Olkhovsky and Recami) dwell time, with 
some numerical evaluations. Then, we pass to examine the equivalence between quantum tun- 
nelling and "photon tunnelling" (evanescent waves propagation), with particular attention to 
tunnellings with Superluminal group-velocities ("Hartman effect"). At last, the main (Cologne, 
Berkeley, Florence, Vienna) experimental data about Superluminal evanescent wave propagation 
are briefly reviewed. 



b : Parte I: TEORIA 



1.0) Introduzione. 

Consideriamo una particella inizialmente libera che, durante il proprio moto, vada ad 
incidere su una barriera di potenziale di energia maggiore della propria. 

Come sappiamo la meccanica quantistica prevede che vi sia una probability non nulla 
che la particella possa attraversare la barriera {effetto tunnel). E dunque lecito chiedersi se sia 
possibile associare un tempo, e quindi una velocita, al processo di attraversamento della barriera 
e, in tal caso, come sia possibile calcolare e misurare sperimentalmente queste quantita. 

Potrebbe sorprendere che una domanda apparentemente cosf semplice e di fondamentale 
importanza nella comprensione di qualunque teoria dello scattering, e non solo, non abbia ancora 

°W Lavoro in parte finanziato da INFN, MURST, CNR (Italia), da CAPES, CNPq (Brasile), e dall' 
I.N.R. (Accademia Ucraina delle Scienze, Kiev). 
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ricevuto una adeguata risposta. 

II problema venne posto per la prima volta nel 1931 da Condon, [1] ed un primo tenta- 
tive di risposta risale all'anno successivo da parte di McColl.[2] Da allora, pero, l'argomento 
rimase quasi ignorato fino agli anni quaranta e cinquanta, con il tentativo di introdurre una 
osservabile quantomeccanica per la variabile tempo, soprattutto nella trattazione delle collisioni 
[per una semplice introduzione in M.Q. di un operatore (non autoaggiunto ma) hermitiano per 
l'osservabile Tempo, si vedano le refs.[3].] 

L'avvento di dispositivi clcttronici ad alta velocita, basati su processi di tunnelling, e 
l'importanza del fenomeno in processi di fissione e fusione nucleare sotto soglia, hanno pero 
suscitato, soprattutto negli ultimi quindici anni, un crescente interesse verso l'argomento, tanto 
da portare, dal 1987 ad oggi, alia pubblicazione di parecchie (almeno una decina[4 — 6]) reviews 
teoriche al riguardo. 

Sperimentalmente, invece, le difficolta nell'effettuare misure non invasive con particelle 
quali ad esempio elettroni, hanno reso molto difficile una verifica dei risultati teorici. Solo negli 
ultimi anni, sfruttando Pequivalenza tra trasmissione di onde elettromagnetiche evanescenti[7] 
e tunnelling di particelle, sono state realizzate alcune misure dei tempi di trasmissione di mi- 
croonde e di fotoni. Ci occuperemo in seguito, quando prenderemo in esame tali esperimenti, 
della suddetta equivalenza. 

Poiche effetti di tunnelling possono intervenire in vari processi fisici (I.teoria dello scatter- 
ing, diseccitazione di stati metastabili, fissione e fusione sottosoglia, etc.) premettiamo che noi 
esamineremo solo il caso unidimensionale di attraversamento, da parte di una particella inizial- 
mente libera, di una barriera di potenziale costante* 1 Vq, che si estenda nell'intervallo [0,d], 
(vedi Fig.I-1); caso peraltro equivalente alle configurazioni degli apparati sperimentali finora 
adottati. 

Malgrado la gran mole di lavori teorici cui accennavamo, non esiste ancora un approccio 
al problema che venga universalmente accettato; noi li raggrupperemo in quattro categorie 
principali. 

Alia prima categoria appartengono tutti quei tempi costruiti "seguendo" il pacchetto in- 
cidente durante l'attraversamento della barriera. 



Figura 1-1. 

Scelta una data caratteristica del pacchetto, per esempio il picco centrale, si con- 
frontano il picco entrante e quello uscente in modo da ricavarne una relazione temporale. Oltre 

0#1 Solo in un caso, quello del tempo di Biittiker-Landauer, la barriera sara anche funzione del tempo. 
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al picco sono stati presi in considerazione dai vari autori anche il centroide (centro di "massa" 
del pacchetto) e il fronte d'onda, molto netto, di un'onda a scalino. 

Tra le critiche piu comunemente avanzate verso questo tipo di approccio, vi e il fatto 
che, secondo molti, il picco uscente non corrisponderebbe sempre a quello entrante della 
presenza, tra le componenti di Fourier del pacchetto, di frequenze corrispondenti ad energie 
prossime o al di sopra deU'energia di barriera. Tali frequenze raggiungerebbero prima delle 
altre la barriera venendo da questa trasmesse piu efficientemente di quelle piii basse. Tutto 
cio, oltre a provocare gli usuali effetti di dispersione, eventualmente presenti anche durante la 
fase di avvicinamento alia barriera, sembra poter avere, in certi casi, un effetto accelerante sul 
pacchetto o, addirittura, portare alia sua trasmissione prima ancora che il picco principale la 
abbia raggiunta. 

A questa categoria di tempi appartiene il phase time. Questo viene ricavato mediante 
l'approssimazione della fase stazionaria utilizzando la definizione di velocita di gruppo del pac- 
chetto e, malgrado le suddette critiche, risulta quello maggiormente in accordo con i dati speri- 
mentali finora a disposizione. 

Un secondo tipo di approccio consiste nell'introduzione di qualche grado di liberta nel 
sistema particella-barriera, per definire un "orologio" che dia una misura del tempo trascorso 
dalla particella all'mterno della barriera. Attraverso l'effetto dell'orologio sulla particella o, 
viceversa, attraverso l'effetto della barriera su un orologio associato alia particella, si risale al 
loro tempo di interazione durante l'attraversamento. 

Biittiker e Landauer, ad esempio, tentano di risalire al tempo di tunnelling calcolando lo 
scambio di quanti di energia da parte di una particella che attraversi una barriera rettangolare 
la cui altezza sia modulata nel tempo (tempo di Biittiker- Landauer) . 

Un altro esempio puo essere la misura della rotazione dello spin di un elettrone quando 
nella porzione di spazio occupata dalla barriera si trovi un campo magnetico uniforme (tempi di 
Larmor) . 

Questo tipo di approccio consente una larga scelta nei gradi di liberta del sistema che pos- 
sono essere usati come orologio e, contemporaneamente, fornisce anche dei metodi sperimentali 
per la verifica dei risultati teorici. Viene pero criticato da alcuni autori, sia perche non tutti i 
tipi di orologio risultano equivalenti, sia soprattutto perche l'introduzione di tali tipi di orologi, 
oltre a modificare il numero di gradi di liberta del sistema, introduce comunque dei processi 
invasivi, i quali condizionerebbero i risultati. Praticamente niente ci assicura che il tempo di 
interazione, in cui comunque lo stato del sistema subisce una perturbazione, corrisponda proprio 
al tempo di attraversamento o di riflessione. 

II terzo tipo di approccio al problema consiste nell' attribuire al moto della particella sotto 
barriera un insieme di traiettorie "semiclassiche" . Con queste viene poi calcolato un tempo medio 
di tunnelling. 
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I cammini possono essere costruiti con vari metodi, per esempio attraverso i path-integrals 
di Feynman, il metodo di Bohm, o l'uso della distribuzione di Wigner: tutti e tre i metodi 
portano, naturalmente ad una distribuzione di tempi. 

Uno degli inconvenienti di questo tipo di approccio e quello di fornire dei risultati comp- 
lessi. Tuttavia e possibile estrarre da questi delle quantita reali (il modulo, la parte reale e la 
parte immaginaria) che, in certi casi, risultano essere strettamente connesse con i tempi definiti 
mediante altri approcci. Proprio per questo motivo l'interpretazione fisica di tali risultati appare 
comunque molto interessante. 

L'ultima categoria di tempi di tunnelling che esamineremo e quella che parte dalla 
definizione di dwell time o tempo di soggiorno. Quest'ultimo viene definite come: 



e rappresenta il tempo speso da una particella albinterno della barriera, o di una qualunque 
regione di spazio, calcolato come rapporto tra densita di probability che la particella si trovi in 
quella regione ed il flusso j- m in essa entrante. 

II problema di tale definizione e che questa ci fornisce si il tempo di soggiorno all'interno 
della barriera, ma senza distinguere tra canale di trasmissione e canale di riflessione. A 
tal proposito, una relazione, spesso considerata ovvia ma non universalmente accettata, che 
dovrebbe legare i tempi corrispondenti ai due canali e: 



Questa relazione, se pure fosse corretta, sembra non bastare da sola ad individuare univocamente 
tt e TR. 

Le critiche a tale relazione, ed i tentativi di alcuni autori di trovare una relazione equiv- 
alente che leghi il tempo trascorso all'interno della barriera con i tempi di riflessione e di 
attraversamento, portano ai cosiddetti approcci spaziali. Questi vengono costruiti attraverso 
l'interpretazione probabilistica standard, in meccanica quantistica, delle densita di corrente rel- 
ative a trasmissione e riflessione. 

Notiamo infine che alcuni dei tempi che descriveremo, potrebbero sembrare dei candidati 
migliori di altri a rappresentare il tempo di tunnelling. Cio sia perche non presentano il com- 
parire delVeffetto Hartman, sia perche in accordo con la (1.0.1). Tuttavia proprio questi tempi 
risultano, invece, maggiormente in disaccordo con i dati sperimentali. 

1.1) Premesse e notazioni. 

Prima di andare a descrivere i phi importanti tipi di tempi ed i relativi risultati premet- 
tiamo alcune considerazioni e notazioni. 




(7.0.0) 



T B = \T(k)\ 2 T T + \R(k)\ 2 T R . 
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Supporremo di aver gia risolto in modo esatto il caso stazionario. Per ogni energia fissata, 
E = h 2 k 2 /2m, sara allora: 

( ijjj = e ikx + R{k)e~ i{kx '^ x < 
^{x-k) = I ^ ll = x {x-k) 0<x<d (1. 1.0) 

{ ip m = T(k)e i( - kx+a ^ x>d 

con R(k) e T(k) ampiezze, rispettivamente, di riflessione e di trasmissione, tali che: 

R(k) = yJl-T(k)*, 

e P = f3(k), a = a(k), corrispondenti ritardi di fase. 
Nel caso particolare di una barriera rettangolare: 



V(x) = { 



V 0<x<d 
altrove, 



x(x;k), R(k), T(k), a(k), j3(k) , sono tutte quantita la cui forma analitica e ben conosciuta; 
avremo dunque: 

f A(k)e~ KX + B(k)e KX E<V 
{ A(k)e- iKX + B(k)e iKX E > V ' 



X(x;k) 



(7.1.1) 



con 



R{k) = < 



T(k) 



. = f y /2m(V - E) /fi, E<V 
" I ^/2m(E - V )/h, E>V 

(k 2 - K 2 ) sinh(Atd) 

[Ak 2 K 2 + (k 2 + K 2 ) 2 sinh 2 (Kii)]2 
(k 2 + k 2 ) sm(nd) 

[4k 2 k 2 + (k 2 - k 2 ) 2 sin 2 (nd)}^ 

2kK 

[4k 2 k 2 + (k 2 + k 2 ) 2 sinh 2 (zed)]* 

2fc« 

I [4k 2 K 2 + (k 2 - K 2 ) 2 Sin 2 (/€<*)] 3 

-2/ck 



a(/e) = < 



arctan 
arctan 

arctan 
arctan 



(fc 2 - K 2 , 

-2kn 

(k 2 + K 2 ) 
'(fc 2 -K 2 ) 

2kK 

\k 2 + K 2 ) 

2kK 



coth(Kd) 



cot(ftd) 

tanh(KeZ) 
tan(ztcZ) 



£ < y 

£ > V 

E < V 
E>V 

E < Vq 
E>V 

E <V 
E>V 



(7.1.2) 



(7.1.3) 



(7.1.4) 



(7.1.5) 



Naturalmente in seguito non avremo a che fare con semplici onde stazionarie ma con 
pacchetti d'onda della forma: 



*(x; t) = J dk Cf(k - k )i>(x; ^e^^ = 
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= J dE g(E - E )^(x; fc(£)) e -*f . (7.1.6) 

Saranno inoltre: 

p=|V(*)| 2 , i = Re {^y,(*)|>(z)}. (7.1.7) 

Nel caso elettromagnetico (microonde e fotoni), che noi supporremo sempre T.E. o T.M., ip 
rappresentera la parte scalare di una delle due conponenti del campo. 

Definiamo infine tempo equivalente, tJ , il tempo che la particella impiegherebbe ad at- 
traversare lo spazio occupato dalla barriera, ma in sua assenza; sara dunque rj q = md/hk. 
Prendiamo ora in rassegna le definizioni di alcuni dei tempi cui abbiamo accennato. 



1.2) Phase Time. 

Supponiamo di avere un pacchetto d'onda molto stretto intorno ad un numero d'onda 
ko. La descrizione della sua evoluzione temporale e in genere molto complicata a causa, spesso, 
della propria natura dispersiva. In ogni caso, sotto opportune condizioni, e possibile seguire la 
posizione del picco di un pacchetto simmetrico con buona precisione, trascurando tali effetti[8]. 
Possiamo quindi provare ad identificare il pacchetto prendendo come riferimento il picco: per 
far cio usiamo il metodo della fase stazionaria. 

II picco del pacchetto sara formato da quelle componenti di Fourier per cui la variazione 
di fase in un intorno di ko sia abbastanza ridotta in modo da non interferire distruttivamente. 
Analogamente il pacchetto trasmesso sara descritto da una funzione d'onda anch'essa dominata 
da una piccola serie di frequenze ciascuna della forma: 

ip(x;k) ~ e i(.kx-Sim+a{k)) _ 

Se vogliamo seguire la posizione del picco, dobbiamo vedere per quali valori di x p (t) la fase e 
stazionaria, cioe, l'argomento dell'esponenziale e massimo.* 2 . Deve quindi essere: 



E(k)t 



, . 1 dE da f i 

Allora a'(ko) = (da/dk)k rappresentera il ritardo spaziale 5x causato dal processo di tunnelling 
e, dividendo per v g (velocita di gruppo del pacchetto d'onda) otteniamo un ritardo temporale: 

Sx ± da fldEX^da da 

6TT = V R =(Vg) dk = {hdk) dk =h d-E- {L2A) 



0#2 Applicando lo stesso ragionamento alia parte incidente otterremmo: x p (t) — j- (^f ) t — (^) t, 
dove (dw/dfc) e proprio, per sua definizione, la velocita di gruppo del pacchetto incidente. Notiamo 
inoltre che, naturalmente, v g = (dw/dfc) dipcndera dalla relazione di dispcrsione u{k) del mezzo in cui il 
pacchetto si propaga. 
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Figura 1-2. 



Notiamo die stiamo calcolando sia v g che tutte le altre quantita in k = ko- Vedremo 
piu avanti se e quando cio possa considerarsi corretto. 

Per motivi che spiegheremo tra poco defmiamo phase time il tempo totale T^(x\,x 2 ;k), 
speso da una particella tra due punti, x\ e X2, esterni alia barriera e abbastanza lontani da essa, 
cioe: ii < e X2 > d. Avremo allora: 

t$(xi,x 2 ; k) = — (x 2 - x x + a'(k)) (1.2.2) 

ed applicando lo stesso tipo di ragionamento alle particelle riflesse: 

T&{x 1 ,x 2 ;k) = -(-2*1 +/?'(&)) (7.2.3) 
v g 

Poiche sia che dipendono linearmente da x\ e x 2 , potremmo pensare di estrapolare 
i tempi di attraversamento e di riflessione direttamente facendo tendere x\ e x 2 rispettivamente 
aOed: cio pero non e corretto. 

Avendo infatti supposto le componenti del pacchetto strettamente distribuite intorno ad 
un numero d'onda ko, se A A; = a, nello spazio ordinario avremo un pacchetto la cui larghezza sara 
dell'ordine di c -1 e, quindi, tanto piu esteso quanto piu esso sara piccato intorno a ko- Dunque 
la parte incidente e quella riflessa della funzione d'onda potranno interferire tra loro anche 
ad una certa distanza (~ o"^ 1 ) dalla barriera (vedi Fig.I-2). Del resto, visto che stiamo usando 
un'approssimazione stazionaria, non stiamo seguendo realmente il pacchetto nella sua evoluzione 
temporale ma, semplicemente, osservando asintoticamente il ritardo di fase corrispondente al 
numero d'onda ko. 

Definiamo comunque ugualmente: 

Ar^ = —{d + a'{k)), 

v g 

Ar £ = -V(fc)), 

Vg 

che chiameremo tempi di fase estrapolati. In seguito, pero, 
carattere puramente asintotico di tali definizioni. 

Nel caso di barriera rettangolare abbiamo: 

Ar%(0,d;k) = A T £(0,d;k) = 



(7.2.4) 
(7.2.5) 

dobbiamo tenere sempre presente il 
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m ( 2ndk 2 (K 2 - k 2 ) + e 4 sinh(2/cd) 
hkn y 4n 2 k 2 + e 4 sinh 2 (nd) 

con e = 2mVo/h. AH'aumentare di d il temine in parentesi tendera a due, e quindi Ar^ e Ar R 
tenderanno a 2m/fikK = 2/v g K: tale valore non dipende affatto dallo spessore della barriera e 
diminuisce all'aumentare di k, cioe per energie molto piccole (k — ► 0). 

Aumentando dunque lo spessore, possiamo trovare delle velocita di tunnelling d/rf ar- 
bitrariamente grandi: effetto Hartman [8] -Fletcher [9] (o semplicemente effetto Hartman). Tale 
effetto, per quanto potrebbe sembrare fisicamente inaccettabile, perche in contrasto con la rela- 
tivita ed il principio di causalita, e stato realmente osservato in tutti gli esperimenti cui abbiamo 
precedentemente accennato. 

Quasi tutti gli autori[10] insistono comunque sul fatto che la violazione del principio di 
causalita e, in questo caso, solo apparente, e che il comparire di velocita Superluminali nei risul- 
tati sperimentali sia in effeti dovuto ad un reshaping* 3 del pacchetto. 



(7.2.6) 



Figura 1-3. 

Didascalia della Fig. 1-3: 

{3a) Coefficienti di trasmissione T(k) calcolati per diversi valori dello spessore d di una barriera 
di potenziale rettangolare di altezza Vq = h 2 e 2 /2m. Nella stessa figura viene inoltre riportato il 
grafico di f(k - k ) = e -(k-k Q ) 2 /2(Ak)^ con k(j = Q ?£) Ak = Q lk(j 

3b) Coefhciente di trasmissione di una barriera di spessore d = 4/e fino ad energie corrispondenti 
a valori di k/e = 5.} 



A questo punto, partendo proprio da un'analisi del phase time, puo essere interessante 
una discussione, puramente qualitativa, degli effetti subiti da tutto l'insieme di frequenze cos- 
tituenti un pacchetto d'onda neH'attraversamento della barriera, per capire se, ed in quali casi, 
possano effettivamente apparire effetti dovuti al reshaping, e se sia possibile evitarli. 

Tale discussione, che non pretende affatto di avere la validita di una dimostrazione, trovera, 
comunque, conferma nei risultati ottenuti seguendo il centroide o l'evoluzione temporale del 
pacchetto, e si rende necessaria perche, come Trp R e Ar^ R , quasi tutte le altre definizioni di 
tempi saranno calcolate in k = /co- 
in Fig.I-3a possiamo vedere rappresentati i valori di T(k), calcolati per diversi valori di d in 
funzione di e. Nello stesso grafico e riportata la funzione di distibuzione di un pacchetto gaussiano 
piccato intorno a ko = 0.7e. f(k — ko) = e~( fc ~ fe °) / 2 ( Afc ) ; con Ak = O.lfeo- Naturalmente il peso 

o#3 p er j-gghaping d e i pacchetto intendiamo il fatto che, attraversando la barriera, le sue componenti 
di energia piii bassa sono, in genere, trasmesse meno efficientemente di quelle di energia piii alta. Cio 
puo provocare un cambiamento della forma del pacchetto nello spazio delle k (una sua accelerazione in 
qucllo delle x), piii o meno evidenti. 
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di ogni componente di Fourier nel pacchetto trasmesso sara dato dal prodotto T(k)f(k — ko). 



Per barriere molto sottili (de <C 1) T{k) e quasi costante e molto vicino ad 1, tranne die 
per valori di k molto piccoli. Purche dunque k /e non sia troppo vicino a 0, la funzione T(k)f(k— 
ko) avra il suo massimo sempre in ko ed il pacchetto uscente non presenters, deformazioni. 

In questo caso il tempo di trasmissione, calcolato sempre come phase time, risulta maggiore 
del tempo equivalente e non ci sono problemi. 



Aumentando lo spessore della barriera, questa comincera a trasmettere sempre "peggio" 
le componenti del pacchetto di energia piu bassa: infatti, per valori di d abbastanza grandi, 
T(k) rimane molto piccolo tranne quando k/e si avvicina molto ad 1 dove, invece, cresce molto 
velocemente. 

Proprio per questo motivo, secondo molti autori, il pacchetto subirebbe un'accelerazione 
in quanto il massimo del pacchetto trasmesso si verrebbe a trovare verso un valore k' > ko- 
Dunque, secondo questi autori, e come se ad attraversare la barriera fossero solo un sottoinsieme 
di frequenze, quelle piu elevate, propagatesi a velocita maggiore di v g , anche prima di raggiungere 
la barriera stessa. In effetti, e difficile che cio avvenga, e puo essere in ogni caso evitato (senza 
che scompaiano le velocita di gruppo Super luminali). Vediamo come. 

Innanzitutto notiamo che la barriera non puo avere alcun effetto amplificante su nessuna 
componente del pacchetto, ma solo fare da filtro per alcune di esse, quindi T(k) e limitata e puo 
salire solo fino ad un valore massimo pari ad uno (come vediamo anche dalle Figg.I-3a e I-3b). 

Per k < e, inoltre, T(k) e strettamente crescente e, nel caso che ci interessa, cioe per 
barriere abbastanza spesse, T(k) cresce molto velocemente solo per valori di k prossimi ad e, 
dove ha un flesso obliquo in cui passa da concava a convessa. Saranno inoltre: 

lim T(k) = . 1 (1.2.7) 

k^e I 



limnfc) = 2^(3 + ^) 
k ^ e 3(l + d 2 e 2 )2 

limiti validi anche per k > e. 
In prossimita di e sara quindi: 

r(fc)a - r J_ + ^ e ( 3 + *? ( fc - e ) (7.2.9) 
7# 3(1 + ^2)5 

e dunque, nel punto in cui cresce piu velocemente, T(k) cresce non piu velocemente di T'(k)(k — 
e). 



Per quanto riguarda f(k — ko) abbiamo invece che: 

f(k - k ) = exp[-(k - k ) 2 /2(Ak) 2 }. 
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L'argomento dell'esponenziale sara minore di uno per (k — ko) < Ak e maggiore di uno se 
(k—ko) > Ak: nel secondo dei due casi f(k—ko) decrescera esponenzialmente con ((k—ko)/Ak) 2 . 



Tranne che per barriere troppo sottili, per le quali abbiamo gia visto che T(k) ~ 1, ci 
aspettiamo che, se (e — ko) ~ Ak, buona parte del tunnelling possa anche essere dovuto alia 
trasmissione delle frequenze di energia piu alta e, quindi, il picco delle k potrebbe effettivamente 
spostarsi in avanti: sara pero, in ogni caso, (k' — ko) ~ Ak. 

Diminuendo, invece, l'energia del pacchetto incidente, prendendo cioe dei ko minori, il 
contributo al tunnelling da parte delle energie piu alte diminuisce sensibilmente in quanto f(k — 
ko) decresce esponenzialmente con {{k — ko)/ Ak) 2 : cio sempre piu velocemente man mano che 
ko si allontana da e. 

Affinche sia possibile un reshaping del pacchetto, tale da provocare uno shifting o il 
formarsi di un secondo picco in avanti, e allora necessario che, in certo intervallo, T(k) cresca 
molto velocemente, abbastanza piu velocemente di quanto decresce f(k — ko) in quello stesso 
intervallo. 

Per l'esattezza dovra essere: 

±[T(k)f(k)] = T'(k)f(k - k ) + T(k)f'(k - k ) > 0, (7.2.10) 

funzione molto difficile da studiare a causa della forma abbastanza complicata di T'{k). 

Poiche pero f'(k — ko) = — f(k - ko), ed f(k — ko) > sempre, la condizione (1.2.10) si 

[Ak) 

riduce a: 

(k - k ) 
(Ak) 



T'(k) - l" 2 o; T(A0 > 0. (I-2.ll) 



Osserviamo allora che, come del resto ci aspettiamo, per k < ko la (1.2.11) e sempre verificata; 
per k > ko invece, affinche il prodotto T(k)f(k — ko) sia crescente, dovra essere: 

T'(k) > T(k)^^. (1.2.12) 

Ma, come abbiamo gia visto, T'(k) e limitata, e per barriere abbastanza spesse ha il suo massimo 
in k = e. Quindi, per quanto possa essere piccola T(k), possiamo sempre trovare dei valori di 
Ak tali che la (1.2.12) non sia piu valida e quindi il picco rimanga in ko- Cio anche per barriere 
molto spesse,* 4 per le quali pero la probabilita di tunnelling diviene infinitesima; basta solo 
che: 1) ko non sia troppo vicino a e, e 2) la distribuzione delle k sia abbastanza stretta intorno 
a ko- 



Figura 1-4. 



°# 4 Per nd = (e 2 - k 2 ) 1/2 d > 1, si ha T(k) = T((e 2 - k 2 ) 1/2 ) Y^^^ ma in tal caso h anche 

T'(k) - e- Kd 
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Didascalia della Fig.I-4: 

{Grafici di T(k)f(k — ko) per energie al di sotto dell'energia di barriera in funzione del numero 
d'onda k, per diversi valori di ko e diversi valori dello spessore a della barriera (ko = 0.3s, 
Fig.I-5a; k = 0.5s, Fig.I-5b; k = 0.7s, Fig.I-5c; k = 0.9s, Fig.I-5d; Ak = 0.1k ; e dall'alto 
verso il basso: a = 1/s, a = 3/e, a = 6/e, a = 10/e, a = 15/e, a = 20/e, a = 25/e. Notiamo 
che i grafici relativi agli ultimi due valori non compaiono nelle prime due figure a causa della 
forte attenuazione del picco.} 



In Fig.I-4 sono riportati i grafici di T(k)f(k — ko), per diversi valori dello spessore 
della barriera e di ko (notare la scala logaritmica) . Come si puo vedere: solo per valori di ko/s 
abbastanza vicini ad uno (ko/s = 0.9) buona parte del tunnel avviene per componenti di energia 
al di sopra dell'energia di barriera. 

Notiamo infine che, anche quando il picco dovesse spostarsi in avanti, cio non avrebbe 
alcuna influenza diretta sulla sua posizione nello spazio delle x, se non fosse che, in tal caso, sia 
Trp R , che Arrp R , perderebbero qualsiasi significato fisico in quanto calcolati in k = ko. 

Tutto cio ci assicura che il pacchetto trasmesso, che andiamo a rivelare, non presenta delle 
grosse deformazioni rispetto a quelle- incidente. 



1.3) Tempi costruiti seguendo il centroide. 



Alcuni autori, anziche seguire il picco attraverso il metodo della fase stazionaria, 
preferiscono usare come riferimento il centroide (o centro di massa) del pacchetto. Cio, sia 
perche, come abbiamo visto, per poter applicare il metodo della fase stazionaria abbiamo bisogno 
di pacchetti molto stretti in k, e questo li rende molto estesi nello spazio ordinario, sia so- 
prattutto perche, cosf facendo, e possibile valutare Peffetto di eventuali accelerazioni causate 
dall'attraversamento della barriera. 



Supponiamo che inizialmente (t < 0) il pachetto si trovi ad una certa distanza dalla 
barriera tale che: 

f oo 

\ib(x,0)\ 2 dx ~ 



che sia, cioe, trascurabile la probabilita che la particella si possa venire a trovare oltre lo 0. 
Identifichiamo la posizione della particella all'istante t = attraverso la posizione del suo "centro 
di massa", cioe: 

, x J°°x\Mx,0)\ 2 dx 

Dato: 

f(k) = f(k,0)-^= J dx^(x,0)e- ik * = \f(k)\e^ (7.3.1) 
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si pub dimostrare che[ll]: 



e ponendo: 



avrcmo: 



x = x(0) = 



fo°°M\f(k)\ 2 % 



< ... >= 



Io°°dk\f(k)\ 2 

/oo 
dk\f(k)\ 2 ; 
-oo 



= - < e'(fc) > 



<...i> 

< ... > in = — — — =< ...1 >; 



< ... > T = 



< ...\T\ 2 > 
< |T| 2 > ' 



< ... >R= 



< 1 > 

< ...\R\ 2 > 

< \R\ 2 > 



Xin(t) = x H < k > in i, (i -> 0); 



m 



= xq H < fe >t t— < a' >t, (i — ► oo); 



£r(*) = xq H < fc > R i- < > R , (t -> oo); 



m 



(1.3.2) 



(1.3.3) 
(1.3.4) 
(7.3.5) 



Facendo evolvere in avanti nel tempo x- m (t) e xr(0), ed indietro xt(£), si possono estrapolare 
iin(0), £r(0) e t"c(d), rispettivamente come tempi per i quali il centroide passa per e per d.# 5 



I tempi di trasmissione e di riflessione saranno dati allora da: 



m 



r£ = t T (d) - tin(0) = - 



r£ = t R (0) - tin(0) 



(eZ-x + <«'>) + 



x 



< k > T 



< fc >i 



-X + < P' >R 



+ 



X 



(1.3.6) 



(1.3.7) 



< k >R < fc >i 

Leavens e Aers[12] dimostrano che, nel limite Ak — > 0, si ha: — > Ar^ e — > At^, e 
che, comunque, le eventuali correzione sono del primo ordine in Ak, il che conferma i ragiona- 
menti di natura qualitativa precedentemente fatti. 



Risultati analoghi sono raggiunti anche da Martin e Landauer[13] che eseguono il calcolo 
nel caso elettromagnetico, e da Collins, Lowe e Barker [3] che seguono l'evoluzione temporale di 
un pacchetto gaussiano attraverso l'equazione di Schrodinger time-dipendent. 

Anche loro perb, pur calcolando esplicitamente il momento in cui il centroide lascia la 
barriera (non c'e autointerferenza per la funzione d'onda per x > d), estrapolano quello in cui 
questa viene raggiunta dal centroide [6]. 

Passiamo ora ad esaminare alcuni dei tempi definiti mediante "orologi". 



°# 5 Notiamo che stiamo sempre usando delle forme asintotiche della funzione d'onda, trascurando 
quindi, sempre, gli effetti di autointerferenza in prossimita della barriera. 
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1.4) Tempi di Biittiker e Landauer. 



Per determinare tt, Biittiker e Landauer[14 — 16], nel 1982, propongono di considerare 
una barriera rettangolare oscillante nel tempo, supponendo il tempo di attraversamento della 
sudetta barriera pari al tempo di interazione delle particelle con il potenziale oscillante. 

Consideriamo allora una barriera di potenziale rettangolare di altezza Vq, a cui venga 
sovrapposto un potenziale oscillante 5V cos cut. 

A frequenze piuttosto basse, il potenziale variera molto lentamente, dunque la particella, 
nell'attraversamento, risentira solo di una parte del ciclo di modulazione e, finche il periodo 
corrispondente all'oscillazione sara lungo rispetto al tempo di attraversamento, la particella 
interagira con un potenziale quasi statico. 

A frequenze piu alte (a; 3> 1/tt), la particella, pur vedendo un potenziale medio Vo, subira 
l'effetto di vari cicli di oscillazione e potra assorbire o cedere dei quanti di energia pari a Tiuj. 

La frequenza a cui avviene la transizione tra il comportamento adiabatico, tipico delle 
basse frequenze, e quello in cui si presenta assorbimento o cessione di energia, fornisce una 
misura del tempo di interazione della particella con la barriera. Naturalmente questa sara 
soltanto un'indicazione approssimata di tale tempo. 



Al primo ordine in 5V appariranno solo le due bande E±%lu. Inoltre le particelle appartenti 
alia banda energetica piu alta saranno favorite nell'attraversamento rispetto a quelle di energia 
piu bassa. 

Biittiker e Landauer mostrano quindi che, per barriere spesse e frequenze non troppo 
elevate (hu; piccolo rispetto sia ad E sia a Vq — E), Pintensita relativa delle due bande sara data 
da: 

I due autori identificano quindi il tempo di attraversamento con: t|? l = md/(hn). 
Nel limite di frequenze molto piccole la (1.4.0) si riduce a: 



T ± {u) 



T{u) 



< 5V 



) • (/.4.1) 



Come ci aspettiamo, il numero di particelle che avranno assorbito o emesso energia e, in 
questo caso, assolutamenete indipendente da w. 
Sempre dalla (1.4.0) si ha : 

= tanh(^ L ) (1.4.2) 

Cio mostra come proprio Tp L determini il passaggio dal comportamento adiabatico a basse 
frequenze, dove T + ~ T_, a quello ad alte frequenze dove T + S> T_. 



Per quanto riguarda le particelle riflesse, sempre nei limiti Tito <C E e %u -C Vq — E, 
trovano: 

R± 2 ,Wt k \1 



R 



V 2% ) 



(1.4.3) 
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con tr = hk/(VoK). 

Notiamo che anche la (1.4.3) e indipendente da oo. 

II comportamento a basse frequenze di eq.(I.4.1) ed eq.(I.4.3) e tipico di un sistema a due 
stati, | 1 > e | 2 >, di energia E ed E±Huj, portati in risonanza da una perturbazione Vicosuit. 
Se per t = l'intera popolazione del sistema si trova nello stato | 1 >, la popolazione dello stato 

V t 2 

I 2 > cresce inizialmente con ( — — ) . Analogamente avviene se le energie dei due livelli E\ ed 

V 2ft/ 

E 2 sono uguali. Nelle (1.4.1) e (1.4.3) r^ L e r^ L giocano quindi effettivamente il ruolo dei tempi 
di interazione del sistema particella-barriera. 



1.5) Tempi di Larmor. 

Nel 1966 Baz'[17, 18] propone di sfruttare la precessione di Larmor, causata dalla presenza 
di un campo magnetico su particelle dotate di spin, per misurare i tempi di collisione di queste 
particelle. 

Nello stesso anno Rybachenko[19] applica questo metodo per calcolare i tempi di tunnelling 
nel caso unidimensionale di una barriera rettangolare. 

Consideriamo allora un fascio di particelle di spin ^, polarizzate in direzione x, massa 
m e energia cinetica E, che si muovono in direzione y (vedi Fig. 1-5). Supponiamo inoltre che 
un campo magnetico omogeneo debole, Bq, rivolto lungo l'asse z, occupi la zona della barriera, 
sovrapponendosi a questa. 



Figura 1-5. 



Seguendo il ragionamento di Rybachenco, le particelle che penetreranno la barriera, 
attraversando il campo magnetico, subiranno una precessione di Larmor con frequenza cjl = 
g/iBo/h, con g rapporto giromagnetico, e fi momento magnetico delle particelle. 

La precessione si arrestera nel momento in cui la particella riemergera da una delle due 
facce della barriera. 
Poiche[20] valgono 

rr ^ I 

< S x > T = - COSW L T , 

^ ■ L 

< by > T = --smuj L T , 
nel limite di campi magnetici deboli avremo: 

n 

< S x > T ~ -, 
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< Sy > T ^ ~^LT^ T , 

dunque, la componente di spin acquistata in direzione y dalla particella sara proporzionale al 
suo tempo di permanenza alPinterno della barriera. Poniamo allora: 

Stranamente, pero, Rybachenko non prende in considerazione l'effetto maggiore del campo 
sulle particelle, cioe l'allineamento degli spin nella propria direzione. Uscendo dalla barriera, 
infatti, gli spin delle particelle avranno acquistato anche una componente lungo z pari a ±h/2. 

Mentre al di fuori della barriera l'energia della particella era indipendente dallo spin, al 
suo interno questa dipendera anche dalla componente z dello spin, a causa dell'effetto Zeeman. 

La differenza di energia per le particelle spin- up e spin-down sara zizhuj-^/2. Le particelle 
spin-up penetreranno quindi meglio la barriera. Ponendo: 

= (l»+| 2 + I-D-I 2 )- 1 ' 2 f° + V<=» 



con: 



D± = T(K ± )e a e~ iK±d 



dove n± e il k corrisponente a E ± — wl. Se: 



< Si > T = - < ij>\a\il> >, 



otteniamo: 



< s, > T = I j^ja + jr^ ' ( / - 5 - 0fl ) 

IT T I 

<S y >T=-ham{a+-a-) + + ~ T 2 , (7.5.06) 



IT T 

< 5a > T = ficos(a + -a_) 12"^ |r 12 • (7.5.0c) 



Espressioni analoghe si trovano per le particelle riflesse, basta sostituire T± con R±. Si 
pub inoltre dimostrare che: 

I J. |2 _ \rp |2 

< S 2 > R = - < S z > T ' + 2 ■ ' ~ ' , (7.5.1a) 

-ft+ + -ft— 



< <5 V >R= — < Sy >T 



i? 1 7?_ 



+ |T_| 



iis + i 2 + i*-i 2 (/ - 5 ' 16) 
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T,T_ 



+ |T_ 



\R+\ 2 + \R\ 



< S x >r= — < S x >T 
Nel caso di campo magnetico debole k± ~ re =F muj^/h, ed inoltre 

| T+ ,2 _ | T i2 TntVL dT 

are dn 



(/.5.1c) 



(1.5.2) 



772 

Nella (1.5.2) il termine - — — , die moltiplica wl, ha, naturalmente, le dimensioni di un 

ILK OK 

tempo. Bxittiker[21] nel 1983 suggerisce allora l'introduzione di tre tempi r^ T , r^ T e r^ T , nel 
seguente modo. Pone: 



< S z > T = (^/2)a; L T^ L T , 



(/.5.3a) 



< S y > T = -{h/2)uj L T^ T , 



(7.5.36) 



<^> T = (V2)[l-(^r^ 2 )/2], 



saranno allora: 



r 2 T = lim n 



< gg >T 

ft, , 
2^L 



lim 



< 5 2 > T 



lim 



< ^ >T 

ft, , 
2^L 



ft, , 

m 

are 



m <91nT 

flK 8k 

m da 
are 8k ' 



(S 2 + 



dlnT\ 2 
8k 



Svolgendo i calcoli Biittiker ottiene: 

me 2 (k 2 - k 2 ) sinh 2 (re<f) + (recfe 2 /2) sinh(2reci) 



are 2 



Ak 2 K 2 + e A sinh 2 (red) 

mA; 2kcI(k 2 — k 2 ) + e 2 sinh(2red) 
are 2 Ak 2 k 2 + e A sinh 2 (red) 



^ L T = V^t + ^t 



Per barriere abbastanza spesse si ha: 



are ' 



t 2mfc 

7- no 

T ^ T ~ ae 2 re 



(/.5.3c) 
(/.5.4a) 
(1.5.46) 
(/.5.4c) 

5.5a 

5.56 
5.5c 

(7.5.6) 



Notiamo dunque che, nel limite di barriere spesse, = r^ L . Del resto, nel caso della 
componente z dello spin, non abbiamo una vera e propria precessione, ma un "salto" in posizione 
spin-up o spin-down, salto accompagnato da uno splittamento dei livelli energetici. 

Lo stesso Biitticher[21] dimostra inoltre che: anche per r^ T si pud fare lo stesso discorso 
del sistema a due livelli fatto per r^ L . Questi due tempi comunque, a mio parere ed in base ai 
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risultati sperimentali, vanno considerati come tempi di risposta della particella alia perturbazione 
in qualche grado di liberta diverso dalla posizione, e quindi non come tempi di attraversmento. 

Per quanto riguarda poi t^ t , e ancora piu difficile assegnargli un significato fisico: infatti, 
se immaginiamo che anche la componente x dello spin preceda intorno all'asse z, allora a tale 
precessione dovrebbe corrispondere una componente media dello spin in direzione x pari a: 

<S X>T = (V2)[l-(c^ L T 2 )/2], 

e non a: 

< S x > T = (ty2)[l - (^r,V)/2]. (/-5.3c) 

Dunque alia luce di quanto detto prima su r^ T , r^ T potrebbe al massimo essere visto 
come una media di r^ T e r^ T e, infatti, alcuni autori introducono direttamente un tempo = 

Nel limite di barriere spesse ~ r^ T . 



L'unico dei tre tempi di Larmor qui definiti, cui sembra quindi possibile attribuire un 
significato fisico, e r^ T . 

Per questo, Falk e Hauge[22] trovano, nel 1988, le seguenti due relazioni: 

t^ t = — (x 2 - xi + a') + ^2 [sin(/3 - 2/cxi) - sin(2a - (3 + 2kx 2 )] (7.5.7a) 

e 

Ty T = —(x 2 -x 1 + a / ) + ^2[sin(/3- 2/cxi) - sin(2a - (3 + 2A;x 2 )] + 

[sin(/3-2/cxi) + sin(2a-/3 + 2/cx 2 )] (7.5.76) 



2hk 2 R [ 

per la parte riflessa , dove x\ e x 2 sono due qualsiasi punti esterni alia barriera (uno alia sua 
destra ed uno alia sua sinistra) e il campo magnetico, anziche essere limitato alia sola zona della 
barriera, e esteso a tutto l'intervallo {x\,x 2 ). 



Si vede facilmente che le (1.5.7) corrispondono ai phase time piii dei termini oscillanti, le 
cui ampiezze aumentano al diminuire dell'energia incidente. 
Nel caso di barriera rettangolare abbiamo: 

T^(d)=T^(d)=T D (d), 

dove r D non e altro che il dwell time, che introdurremo in seguito. 



1.6) Tempi complessi: path-integral. 

L'introduzione di tempi complessi nasce inizialmente dalla seguente idea: al di sopra 
dell'energia di barriera v = Hk = hy/ k 2 — e 2 e, dunque, tt = - = -z — • 

V JlK 
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Per E < Vq, invece, il vettore d'onda diviene immaginario; consideriamo pero lo stesso il 
moto della particella sotto barriera, come se esso avvenisse lungo una traiettoria classica, ma a 
velocita e tempo immaginari. 

Naturalmente non e possibile dare a queste quantita alcun significato fisico, in quanto, 
se "dall'entrata" della particella all'interno della barriera, alia sua "uscita" da una delle due 
faccie, sara trascorso un certo tempo, per quanto piccolo (o grande) questo possa essere, dovra 
per forza essere una quantita reale. 

Potremmo allora immaginare che siano le traiettorie ad essere complesse, ed i tempi e le 
velocita reali, considerando le quantita: 



H\k\ 



m 



s d md 
Tt v Tin' 

dove la S sta per semiclassiche. 

Naturalmente anche e fisicamente inaccettabile, in quanto diverge per k = e e 
resterebbe poi, in ogni caso, il problema delle particelle riflesse. 

Notiamo comunque che, nel caso di barriere spesse anche r t BL , e naturalmente r^ T , tendono 

Per quanto fisicamente inaccettabile, e interessante notare che nel 1992 Hagmann propone 
di considerare il caso di una particella che, per attraversare la barriera, acquisti una certa 
energia AE, per un intervallo di tempo At. II risultato che viene fuori applicando il principio 

.. . . . , . . c md 

di mdetermmazione e proprio At = t£ = — — . 

Tin 

Tornando ai tempi complessi, nel 1987, Sokolovski e Baskin[23] propongono una general- 
izzazione del concetto classico di tempo alia meccanica quantistica ed applicano, poi, il metodo 
proprio al caso del tunnelling. 

Cosideriamo una particella che, emessa all'istante t\ nel punto ri, sia rivelata all'istante 
ti nel punto f 2 . Supponiamo inoltre che la particella, muovendosi lungo la traiettoria r(t), in 
un potenziale V(f), abbia attraversato una certa zona di spazio Q. 

II tempo trascorso dalla particella in quella regione di spazio sara dato allora da: 

tS = f 2 dte n (m), (J-6.0) 
Jti 

dove @n(f(t)) e 1 se f(t) appartiene ad Q, altrimenti. Nel caso unidimensionale avremo: 

= f 2 dt [ d dx 5(x - x(t)). (7.6.1) 
Jti Jo 

Se allora usiamo il metodo del path-integral di Feynman per costrurci delle traiettorie su 
cui mediare i tempi, otteniamo: 

T n (xi,t 1 ;x 2 ,t 2 ;k) =<t2[x(.)] > pat h, 
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Re r" = Ty T , (/.6.3a) 



in cui x(.) e un cammino arbitrario (nello spazio delle fasi) tra (xi,t{) e {0:2,^2) ■ In generale r n 
sara complesso. 

Sokolovski e Baskin trovano quindi: 

con A = Te ta , B = Re 1 ^ . Nello stesso articolo i due autori trovano una relazione tra e i 
tempi di Larmor. 

r 2/T) 

Im = r^ L T , (7.6.36) 

|r£| = -r x L T . (/.6.3c) 
Analogamente per la parte riflessa avremo: 

Re r£ = T^p; , (/.6.4a) 

Im r£ = r z L R , (Z.6.46) 
|t^|=t x l r . (/.6.4c) 

Malgrado la connessione cosi sorprendente tra r n e i tempi di Larmor, e difficile dare 
una interpretazione fisica di tali risultati, in quanto, come gia detto, se esiste un tempo legato 
all'attraversamento della barriera, questo deve comunque essere reale. 

Una (apparentemente un po' assurda, ma forse non tanto quanto sembra: si ricordi ad es. 
la teoria dell' elettr one, con cronone, di Caldirola: cfr. R. H.A.Farias e E.Recami, "Introduction of 
a quantum of Time ('chronon') and its consequences for quantum mechanics" , archivi elettronici 
LANL # quant-ph /9706059| ) e forse quella data da Hanggi[24] nel 1993. Secondo quest'autore 



il tempo di tunnelling sarebbe caratterizzato da due scale di tempi. Non sa pero spiegare dare 
un significato fisico all'idea. 

Sokolovski e Connor, [25] nello stesso anno, criticando tale affermazione, insistono sul fatto 
che e impossibile e nsicamente insensato pensare che vi possano essere due tempi di attraversa- 
mento diversi. Concludono quindi che quello che va considerato come tempo di attraversamento 
e il modulo di . 

Soffermiamoci, pero, un attimo sull'ipotesi di Hanggi. Se guardiamo alia forma dell'onda 
trasmessa, possiamo osservare che: 

= ^j m ( x , k) = T(k)e ia e ikx = e^^eV. (1.6.5) 

Notiamo allora che ogni qual volta applichiamo certi metodi otteniamo un tempo legato 
a da/dE, ed uno legato a dhiT/dE. 
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Potremmo allora immaginare questo secondo tempo come il tempo che impiega la barriere 
a smorzare il segnale durante l'attraversamento, tempo che non influirebbe affatto sulla sua ve- 
locity. Appare comunque strana, in questo caso, la dipendenza di tale tempo da d, in quanto, 
in questo modo, anziche diminuire all'aumentare dello spessore della barriera, aumenta pro- 
porzionalmente a esso. 

Andrebbe poi spiegato come tale tempo possa essere legato alle transizioni di livello, sia 
nel caso di Biittiker-Landauer, sia nel caso degli spin. 

1.7) Tempi complessi: metodo di Bohm. 

Sempre nell'ambito dei tempi complessi, Leavens e Aers[26], nel 1993, partendo dallo 
stesso operatore introdotto da Sokolovski e Baskin (r^), propongono l'uso del metodo di Bohm 
per ricavare delle traiettorie "semiclassiche" da adoperare per il calcolo del tempo medio. 

II metodo di Bohm, come sappiamo, ci fornisce un'equazione del tutto equivalente 
all'equazione di Schrodinger, consentendoci contemporaneamente un'interpretazione phi "clas- 
sica" dei fenomeni quantistici. Accenniamo brevemente ad esso. 

Poniamo V = Re tS ^ h , con R, S reali. Applicando l'equazione di Schrodinger a ip e sep- 
arando la parte reale e quella immaginaria, otteniamo una equazione di continuita (la seconda 
delle due), ed un'equazione equivalente a quella di Hamilton- Jacobi per un potenziale V(x) 
modificato, a cui e stato aggiunto il terminer 



e la cui soluzione ci da S. 

La R invece viene ricavata dall'equazione di continuita. 

Se consideriamo l'approssimazione WKB, e trascuriamo il potenziale quantistico Q, otteni- 
amo per la S delle soluzioni di un'equazione di Hamilton-Jacobi classica nel potenziale originario: 
naturalmente, in questo caso, S pub essere, anzi sara spesso, una quantita complessa. Notiamo 
infine che, anche se quelle che otteniamo con il metodo di Bohm son traiettorie classiche, queste 
provengono da un potenziale modificato. 

I risultati fin qui ottenuti con questo metodo portano comunque a risultati fisicamente 
assurdi (vedi ref. [6] e [24]) e quindi non ce ne occupiamo. 

1.8) Dwell time. 

II Dwell time viene introdotto per la prima volta da Smith[27] nel 1960, alio scopo di 
calcolare la durata media di un processo di collisione senza distinguere tra i vari canali e, come 
gia detto, viene definite come rapporto tra la probability che la particella si trovi in una certa 
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regione di spazio ed il flusso j- m entrante in quella stessa regione, senza considerare, nel nostro 
caso, se la particella verra riflessa o trasmessa: 

T B (xi,x 2 ;k) = jr 1 / \ifj(x,k)\ 2 dx = 

J X\ 

fX2 

/ \ip(x,k)\ 2 dx. (7.0.0) 

J X\ 



V 

Nel caso di barriera rettangolare abbiamo: 

.2 _ r,2\ , c 2 



n . , mk 2Kd(K 2 — k 2 ) + e 2 smh(2Kd) 

r D (xi,x 2 ;k) = — r — '— ^ — i - (7.8.0 

v ; hn 2 4fe 2 K 2 + e 4 sinh 2 (Kd) V 7 



che, per kc(> 1, diventa: 

r B = = ^ (/.8.1) 
Vqk he z K 



Per barriere abbastanza spesse, dunque, anche r d , come r^ R e Arrp R , diventa indipen- 
dente dal loro spessore ma, in forte contrasto con questi, diminuisce al diminuire di k, fino ad 
annullarsi per k = 0. 

Data la sua definizione, ed il fatto che la trasmissione o la riflessione di una particella da 
parte di una barriera sono eventi mutualmente esclusivi, molti autori (quasi tutti) concludono 
che qualunque tempo candidato a rappresentare tt e tr debba necessariamente verificare la 
relazione: 

Tv = \T(k)\ 2 T T + \R(k)\ 2 T R . (7.0.1) 

Tale relazione e ad esempio verificata da r n . Sokolovski e Baskin[23] nel 1987 trovano 
infatti la relazione: 

r o = + | T |2 r n_ 

Da questa separando parte reale e parte immaginaria si ottiene poi: 

i^ = \R\ 2 r^ + \T\ 2 r^ 
\\R\ 2 t^+\T\ 2 t^ = 0. 

La prima di queste due equazioni viene ricavata indipendentemente anche da Falck e 
Hauge[22] nell'anno successive La seconda, invece, non ci da altro che una legge di conservazione 
del momento angolare, o meglio, in questo caso, della componente z dello spin. 

La (1. 01) invece sembra non essere verificata nel caso dei phase time e, a maggior ragione, 
dei tempi di fase "estrapolati" . Anzi, Hauge et al.fll] trovano la seguente relazione: 

T D (xi,x 2 ;k) = 

\T(k)\ 2 r%(x 1 ,x 2 - k) + \R(k)\ 2 r^ Xl ,x 2 ; k)+ 

^ S m([3-2kx 1 ) (7.8.2) 
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Tuttavia, i due autori dimostrano anche che se consideriamo che ogni pacchetto ha una 
certa larghezza in k (Ak = a) e applichiamo la (1.8.2) a tutto il pacchetto, abbiamo:* 6 

< t b (xi, X 2 \ k) >~ 

< \T(k)\ 2 >< T%( Xl ,x 2 ;k) > + < |i?(£;)| 2 >< T*{ Xl ,x 2 ;k) > + 

i^ff- 1 / dk sin(/3 - 2kxi) + 0(a) (7.8.3) 

Se dunque \x\\ ^> cr _1 , l'argomento dell'integrale oscillera abbastanza velocemente da renderlo 
trascurabile e quindi possiamo scrivere: 

< t d (xi,x 2 ;£0 >= 

< \T(k)\ 2 >< T%( Xl ,x 2 ;k) > + < \R(k)\ 2 ><T&(x 1 ,x 2 ;k) > (7.8.4) 
Rclazione in accordo con la (1.0.1) a meno di 0(a). 

Naturalmente la (1.8.4) non e valida per i tempi di fase estrapolati, e mette nuovamente 
in evidenza il carattere puramente asintotico dei phase time. 

Prima di passare alle critiche mosse da alcuni autori nei confronti del dwell time, analizzi- 
amo un po' meglio la (1.8.2). 

Secondo Hauge e Stovneng[5] questa proverebbe che r D rappresenta proprio il tempo 

esatto trascorso dalle particelle alPinterno della barriera, mentre il termine — —r sin(/3 — 2kx±) 

nk z 

dovrebbe rappresentare un At provocato dagli effetti di auto-interferenza. 

Infatti, se calcoliamo il dwell time in un'intervallo (— L, x\) e facciamo tendere 7 
all'infinito, il valore che otteniamo divergera anch'esso a +oo. Sottraendo pero il dwell time cal- 
colato in (— L,x\) solo sulla particella incidente, otteniamo* 7 un Ar D (x < Xi;k) che, facendo 
i conti, risulta proprio: 

At d (x < Xl ;k) = sin(/3 - 2 Xl ). (7.8.5) 

La (1.8.2) puo quindi essere riscritta come: 

\T(k)\ 2 4( Xl ,x 2] k) + \R(k)\ 2 Tg(x 1 ,x 2 ;k) = 
t d (xi,x 2 ; k) + Ar D (x < x±; k) = 

T D (xi,x 2 ;k)-^sm(P-2kx 1 ) (7.8.6) 



°# 6 Naturalmente in (1.8.3) dovremmo considerare < \T(k)\ 2 r^ (xi, X2\ k) > e 
< \R(k)\ 2 T^{x 1 ,x 2 -k) >, invece di < |T(fc)| 2 >< t$(x u x 2 ; k) > e < \R(k)\ 2 >< T%(x 1 ,x 2 ;k) >. Si 
puo pero dimostrare[4] che l'errore che si commette usando la (1.8.4) e dell'ordine di a. 

Un procedimento simile, in cui si valutano pero i flussi positivi e quelli negativi separatamente, 
verra adottato da Olkhovsky e Recami per generalizzare la definizione di dwell time. 
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A favore di questa tesi i due autori portano il fatto che il suddetto termine di auto- 
interferenza e del tutto indipendente da T(k) e da a(k), proprio perche non c'e interferenza per 
x > d. Aggiungono quindi gli esempi di due casi limite. 

II primo riguarda una barriera infinitamente spessa (in pratica uno scalino). Essendo la 
barriera infinitamente spessa, non ci saranno particelle trasmesse e tutte le particelle verranno 
riflesse: R = 1. E allora facile dimostrare che: 

M = - = |? 

kv nkn 

In questo caso scompare l'apparente contraddizione tra il tempo di fase estrapolato, che au- 
menta con k~ l al tendere di k — ► 0, e il dwell time che invece va a con E/v g ~ k. Infatti, se r D e 
il tempo speso all'interno della barriera e Ar D il ritardo (o anticipo) dovuto all'autointerferenza, 
poiche al diminuire dell'energia incidente la particella penetrera sempre meno all'interno della 
barriera, a prevalere sara il secondo termine, cioe Ar D . 

II secondo esempio riguarda, invece, una barriera a forma di delta di Dirac. Questo caso 
e uno dei primi ad essere trattati e risolti nella letteratura sull'argomento. Naturalmente, in 
questo caso, il dwell time sara nullo, ma si ha:* 8 

At£ = At£ = T(*)-W 

con: 

T{k) = 777? - 

II tempo di tunnelling non puo quindi che provenire dal termine di autointerferenza. 



1.9) Generalizzazione del dwell time. 



Come gia accennato non tutti gli autori concordano sull'importanza fin qui assegnata 
al dwell time anche sulla base, ma non solo, della (1.0.1). Tale relazione, infatti, invocata 
come conseguenza del principio di sovrapposizione, e della mutua esclusivita degli eventi di 
trasmissione o di riflessione, implicherebbe: 

/ |V(x,fc)| 2 dx = j(|TW| 2 r T + |#)| 2 r R ) = 

J Barriera 

= \T(K)\ 2 jr T + \R(k)\ 2 jr R (7.9.0) 

o#8 j n q ues to caso d ha solo la funzione di un parametro, in quanto la S di Dirac e costruita come 
limite di barriere sempre phi strette ed alte ma di area costante Vod. 
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die quantomeno impone, a priori, tt = tr, indipendentemente dalla forma della barriera di 
potenziale. 



Oltre a cio, la (1. 0.0) viene ricavata da Biittiker[21] nel 1983 seguendo un metodo crit- 
icato, invece, da Olkhovsky e Recami.[6, 28, 29] Infatti pur partendo da un'espressione in cui 
appaiono esplicitamente le ip(x, t), nella definizione non si terrebbe realmente in conto della reale 
evoluzione temporale del pacchetto. Inoltre, a parte la relazione (1.0.1) la suddetta definizione 
non suggerisce nessun metodo per la separazione dei tempi dei diversi processi. Per questo mo- 
tivo, nel tentivo di separare questi tempi, Olkhovsky e Recami introducono preliminarmente le 
seguenti definizioni relative ai tempi di trasmissione e di riffessione: 

-in 



T T = t(xf) T - t(xi) in = 

JZodtJ^t) JZodtJUxut) ~ 

S^dEv\g(E)T\ 2 TP h (xi,x f ;E) _ 
Jg°dEv\g{E)T\* 

= (x f - Xi) < v' 1 > T + < Stt > t (/.9.0a) 



TR = t(x t ) R ~ t{ Xi ) in = 

I-oodttJ l R (x { ,t) JZodttJ^xut) 

' I-oodtJU^t) f-oodtJin(x U t) ~ 

J^dEv\g(E)R\ 2 rP h ^i,xf,E) 
J~dEv\g(E)R\* 

= (x f -Xi)< v- 1 > R + < Str > r (7.9.06) 

Infatti, poiche J(x,t)dt rappresenta la densita di probabilita che una particella passi per 
il punto x, nell'intervallo di tempo (t, t + dt), per determinare il tempo medio in cui un pacchetto 
d'onda ^f(x,t) passa per il punto x, dobbiamo fare la media della variabile t pesata attraverso: 

w(x,t)= f00 J( .f| w . (7.9.1) 

Subito dopo, pero, gli stessi autori[6] notano che le definizioni (1.9.0) sono valide solo 
quando i pacchetti d'onda incidente e trasmesso sono completamente separati sia spazialmente 
che temporalmente. Infatti quando x\ ed Xf non sono abbastanza lontani dagli estremi della 
barriera e possibile avere effetti di interfernza tra la parte incidente e quella riflessa. 

Inoltre, la densita di corrente J(x, t) puo, in genere, cambiare segno durante l'evoluzione 
temporale del pacchetto (per esempio quando il picco dell'onda incidente raggiunge la barriera), 
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cosicche gli integrali J^° OQ dttJ(x,t), che rappresentano la somma algebrica di quantita (flussi) 
positive e negative, e le densita di probabilita w(x,t), potranno non essere phi definiti positivi. 

In tal caso, ciascuna densita di probabilita acquista un significato fisico immediato solo 
durante gli intervalli di tempo in cui la corrispondente densita di corrente non cambia direzione: 
occorre dunque spezzare Pintegrale precedente in vari integrali, ciascuno dei quali sia considerato 
su un intervallo di tempo in cui il segno di J(x, t) sia solo positivo o solo negative In tal modo 
si otterranno delle densita di probabilita tutte definite positive: 

J + (x,t)dt 

w + (x,t) - 



u>_(x, t) 



J_(x,t)dt 



dove J + e J_ rappresentano, rispettivamente, i valori positivi e negativi di J(x,t). 

Alia luce di queste osservazioni, i due autori[6] propongono come tempi medi di trasmis- 
sione e di riflessione le seguenti espressioni: 



TT = t{X{) + - t(Xi) + = 

!^dtJ + {x u t) JZodtJ+ixut) 



(J.9.2o) 



TR = t(Xi)_ - t(xi) + = 
J^dttJ^t) J^dttJ+jxut) 

sZodtj-izut) r QO dtj + (x i ,t) 

Prima di passare oltre, partendo dall'equazione di continuita: 



(7.9.26) 



dp(x,t) dJ(x,t) = Q 
dt dx 

e attraverso l'interpretazione probabilistica standard della p(x,t), vogliamo dimostrare come le 
w±(x, t) appena definite corrispondano proprio alia probabilita che la nostra particella, muoven- 
dosi in avanti o venendo indietro, passi nell'intervallo di tempo (t, t + dt) per il punto x. 

In ogni intervallo di tempo in cui J = J + o J = J_, possiamo scrivere l'equazione di 
continuita applicandola a J± [l'equazione di continuita e in ogni caso valida sempre]: 

d P> (x,t) _ dJ+(x,t) (I 9 3a) 

dt dx 

d P < {x,t) = _dJ (s,t) 
at ox 
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ottenendo cosf le due quantita, dp > (x,t)/dt e dp > (x,t)/dt. Integrando quindi rispetto al tempo 
in un intervallo (— oo,i) possiamo definire: 

P>(x,t) = - fJJ±^dt' (/.9.4a) 

p< ( x ^) = - fj-^^dt' (/.9.4a) 

Imponiamo anche la condizione che p > (x, — oo) = 0, p K (x, — oo) = 0: supponiamo cioe che inizil- 
mente la particella (o il pacchetto d'onda) sia infinitamente lontano da x. Integrando nuovamente 
rispetto ad x otteniamo altre due quantita che indicheremo con N > (x, oo; t), N K (— oo, x; t) e per 
le quali si ha: 

[•CO 

N > (x,cc;t) = / p > (x',t)dx' = 

Jx 



L 



t 

J+(x,t')dt' > 0, (/.9.5a) 



N < (—cc,x;t)= / p < (x',t)dx' 



— CO 



f 

J — . 



J_(x,t')dt' > 0. (7.9.56) 

-oo 

Queste ultime due espressioni ci daranno la probabilita che la nostra particella, muovendosi 
in avanti o all'mdietro, si trovi, al tempo t, rispettivamente a destra o a sinistra del punto x in 
funzione delle densita di corrente J±(x,t). 

Notiamo che le condizioni p > (x, — oo) = 0, p < (x, — oo) = che avevamo imposto prima di 
integrare equivalgono ora a J±(—oo,t) = 0. 

Finalmente, differenziando nuovamente le (1.9.5), questa volta rispetto al tempo, abbiamo: 

d 

J+(x,t) = -j^N > (x,oo;t),>0 (/.9.6a) 
J_(x,t) = —N < (-oo,x;t),> (7.9.66) 

e dunque: 

d 

—N^x^it) 

w + (x,t) = -g-. r (/.9.7a) 

+ v ' ^(^,-00,00) v ' 

d 

—N^-oo^x-t) 

w-(x,t) = ^— -. (/.9.7a) 

v ; ^(-00,2;, 00) v ; 

Tali relazioni bastano a giustificare la suddetta interpretazione probabilistica delle w + (x, t) 
e W-(x, t). 



A questo punto possiamo definire il valor medio dell'istante in cui la particella passa per 
il punto x muovendosi nella direzione positiva o negativa delPasse. Saranno dunque: 
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f°° t J+(x,t)dt 

/" tJ_(x,t)di 



'-<*> - E^-(x,t)dt ( "- 86) 

Abbiamo ormai inoltre tutti i mezzi per poter definire anche le varianze delle distribuzioni 
relative ai suddetti tempi, e saranno: 

Siamo quindi riusciti a costruire un formalismo che ci consente di ricavare sia i valori 
medi, che le varianze (o eventuali altri momenti), relativi alle "distribuzioni temporali" di tutti i 
possibili processi relativi al tunnelling, nel caso unidimensionale. Le stesse definizioni, comunque, 
sono estendibili a qualunque altro processo di collisione anche diverso dal tunnelling e con 
qualsiasi tipo di potenziale. 

Come abbiamo gia visto, per i tempi di tunnelling e di riflessione abbiamo: 



T T (xi, x f ) = t(x { ) + - t(xi) + = (/.9.2a) 
con — oo < x\ < e d < Xf < oo e, in base alle (1.9.9): 

a 2 (r T (xi,x f )) = a 3 (t + (x t )) + a 2 {t + (x i )). (7.9.10) 

Nel caso inolte x\ = 0, X{ = d, abbiamo: 



(v^(0,d)=t(d) + -t(0) + , . 

\a 2 (r T (0,d))=a 2 (t + (d))+a 2 (t + (0)). [ ' ' > 

Prendendo per esempio x\ = e < Xf < d, possiamo ricavare tempi di penetramento 
all'interno della barriera come: 



t^(0, x f ) = t(x { ) + - t(0) + (7.9.12) 

o, per < x < d: 

tr^~(x, x) = t(x)_ - t(x) + (7.9.13) 

e, per — oo < x\ < d: 

TR{X uXi ) = t(xO + - t(X0 + . (7.9.14) 



Infine riesaminiamo, in base alia definizioni qui riportate, quelle precedentemente date di 
phase time e dwell time. 
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Per quanto riguarda il primo, salta per Pennesima volta fuori il suo carattere puramente 
asintotico. Infatti, essendo questo ricavato in un contesto esplicitamente stazionario Tunica 
situazione in cui gli si pub dare un senso fisico in base ai risultati qui riportati e quando x\ — ► 
oo, quando cioe, J + (x,t) e la densita di corrente del pacchetto iniziale in assoluta assenza di 
interferenza, tra la parte trasmessa e quella riflessa, dovuta alia onde riflesse. Analogamente, il 
dwell time, rappresentato dall'espressione equivalente:[30, 31] 



r D (xi,x f ) 



f 

J — i 



i -i 



J- m (xi,t) dt 



t J(x t ,t) dt- t J(x u t) dt 

-oo J — oo 

con — oo < x; < 0, and Xf > d, non e in genere fisicamente significativo: infatti il peso nelle 
medie temporali e definito positivo, e normalizzato ad 1, solamente nei rari casi in cui x\ — ► — oo 
and Ji n = Jin (i.e., quando la barriera e trasparente) . 



1. 10) Tempi di penetrazione e di ritorno: risultati numerici. 

Le (1.9.2) e le (1.9.11-14) non consentono un calcolo analitico agevole per ricavare delle 
espresioni per i tempi di tunnelling (penetrazione, ritorno e riflessione) , nemmeno nel caso 
alquanto semplice di barriera rettangolare. 

Presentiamo dunque i risultati di alcuni calcoli numerici relativi alia durata media di vari 
processi di penetrazione e ritorno di pacchetti gaussiani all'interno di una barriera rettango- 
lare, svolti da Olkhovsky et al. in un recente articolo[29]. Tali calcoli, comunque, confermano 
l'apparire dell' effetto Hartman e, in base a verifiche indirette, risultano, secondo gli autori, ab- 
bastanza in accordo con i dati sperimentali di Colonia, Berkeley, Firenze, Vienna. 



Figura 1-6. 

Didascalia della Fig.I-6: 

{6a) Andamento del tempo di penetrazione medio, rp cn (0, x) (espresso in secondi), in funzione 
della lunghezza di penetrazione Xf = x (espressa in in Angstrom) per una barriera rettangolare 
di ampiezza d = 5 A, con Ak = 0.02 A (linea tratteggiata) e Ak = 0.01 A (linea con- 
tinua). Vale la pena di notare che Tp^(0, x) cresce rapidamente per i primi Angstrom (~ 2.5 
A) iniziali, e tende poi ad un valore di saturazione. Questo sembra confermare l'esistenza del 
cosiddetto "effetto Hartman". 

6b) Come in 6a) con Ak = 0.01 A , ma con una barriera di spessore doppio d = 10 A. Os- 
serviamo che i valori numerici del tempo di tunnelling totale rT un (0, d) restano praticamente 
inalterati quando si passa da d = 5 A a d = 10 A, ad ulteriore conferma delkapparire 
dell'effetto Hartman.} 
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Ricodiamo che, in base alle notazioni precedentemente riportate, e 

*inOM) = / Cf(k - k) exp[ikx - iEt/fh] dk 
Jo 

dove: 

(fc-fc) 2 

f(k-k) = e W 1 , 

E = h 2 k 2 /2m, C e una costante di normalizzazione ed m e in questo caso la massa dell'elettrone. 
Le lunghezze di penetrazione saranno espresse in Angstroms, e il tempo di penetrazione in sec- 
ondi. 

In Fig.I-6a sono riportati i grafici di rp cn (0, x), < x < d, corrispondenti a d = 5 A, 
per Ak = 0.02 e 0.01 A . Come si puo vedere il tempo di penetrazione rp en (0,x) tende 
sempre ad un valore di massimo di saturazione. 

In Fig.6b mostriamo invece il grafico corrispondente a d = 10 A, Ak = 0.01 A . E inter- 
essante notare che, a parita di Ak, rp en e praticamente lo stesso, sia per d = 5 che per d = 10 A, 
un risultato questo che conferma ancora una volta l'apparire del cosiddetto effetto Hartman.. 
Risultati analoghi sono stati ottenuti anche per d > 10 A, facendo variare il parametro Ak tra 
0.005 e le energie E nel range da 1 a 10 eV. 

Figura 1-7. 

Didascalia della Fig.I-7: 

{Andamento di rp en (0, x) (espresso in secondi) in funzione di x (in Angstroms), relativo al tun- 
nelling attraverso una barriera di spessore d = 5 A e per vari valori di E e di Ak: curva 1: 
Ak = 0.02A" 1 e £ = 2.5 eV; curva 2: Ak = 0.02A" 1 eE = 5.0 eV; curva 3: Ak = 0.02A" 1 
eE = 7.5 eV; curva 4: A A; = 0.04A" 1 el = 5.0 eV.} 

Nelle Figg.I-6, 1-7, 1-8 sono riportati gli andamenti delle durate medie per i processi di 
penetrazione e ritorno in funzione della lunghezza di penetrazione (con Xi = 0e0< xj = x < d), 
per barriere di altezza Vq = 10 eV e ampiezza d = 5 A o in alcuni casi 10 A. In particolare: 

- In Fig. 1-6, sono riportati i grafici di rp en (0, x) correspondent! a diversi valori dell'energia 
cinetica media: E = 2.5 eV, 5 eV e 7.5 eV, Ak = 0.02A (curve 1, 2 e 3); E = 5 eV, 
Ak = 0.04A -1 (curva 4); 

per tutti e quattro i casi e d = 5A. 

- In Fig.I-7 mostriamo i grafici di rp CI1 (0, x) corrispondenti a: d = 5A, con Ak = 0.02A e 
0.04 A" 1 (curve 1 and 2); d = 10A, con Ak = 0.02 AA^ 1 e 0.04 A -1 (curve 3 e 4,); dove 
l'enegia cinetica media E e 5 eV, e cioe meta dell'energia di barriera Vq. 

- Infine in Fig.I-8 mostriamo alcuni grafici di TR c t(x, x). Le curve 1, 2 e 3 corrispondono rispet- 
tivamente a: E = 2.5 eV, 5 eV e 7.5 eV, per Ak = 0.02A" 1 e d = 5A; 

le curve 4, 5 and 6 corrispondono invece a: E = 2.5 eV, 5 eV e 7.5 eV, per A A; = 0.04A e 
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d = 5A; 

mentre le curve 7, 8 e 9 corrispondono rispettivamente a: Ak = 0.02A e 0.04 A , E = 5 
eV, d = lOA. 

Figura 1-8. 

Didascalia della Fig. 1-8: 

{Andamento di rp cn (0, x) (in secondi) in funzione di x (in Angstroms) per E = 5 eV, e per vari 
valor i di d e Ak: curva 1: d = 5 A, Ak = 0.02A ; curva 2: d = 5 A, Ak = 0.04A ; curva 
3: d= 10 A, Afe = 0.02A _1 ; curva 4: d = 10 A, Ak = 0.04A -1 .} 

Riguardo al modo in cui sono stati effettuati questi calcoli, gli autori fanno osservare che 
l'integrazione su di e stata eseguita usando l'intervallo temporale [— 10 -13 , +10~ 13 ] s, simmetrico 
rispetto a t = 0. Tale intervallo risulta di tre ordini di grandezza maggiore dell'estensione tem- 
porale del pacchetto d'onde che ricordiamo essere dell'ordine di l/(vAk) = (Ak ^JlE/m)^ 1 ~ 
10~ 16 s. 

Sottolineamo che cio equivale a considerare l'evoluzione del pacchetto d'onda attraverso 
un intervallo temporale [—00, +00], e non invece come se questo iniziasse a evolversi ad un 
certo istante t finite Tutto cio in accordo con le relazioni J±(— 00, t) = o equivalentemente 
p K (x, —00) = 0. Aggiungiamo inoltre che il pacchetto e costruito in modo tale da arrivare con 
il centroide in xo al tempo t = 0. 

Figura 1-9. 

Didascalia della Fig.I-9: 

{Andamento di tr^(x, x) (in secondi) in funzione di x (in Angstrom) per diversi valori di d, E 
e Ak: curva 1: d = 5 A, E = 2.5 eV e Ak = 0.02A" 1 ; curva 2: d = 5 A, E = 5.0 eV 
e Ak = 0.02A -1 ; curva 3: d = 5 A, E = 7.5 eV e Ak = 0.02A" 1 ; curva 4: d = 5 A, 
E = 2.5 eV e Ak = 0.02A" 1 ; curva 5: d = 5 A, E = 5.0 eV e Ak = 0.02A" 1 ; curva 6: 
d = 5 A, E = 7.5 eV e Ak = 0.02A" 1 ; curva 7: d = 5 A, E = 5.0 eV e Ak = 0.02A" 1 ; 
curva 8: d = 5 A, E = 5.0 eV e Ak = 0.02A~\} 

Dalle suddette Figure 1-7)— 1-9) si pud vedere che: 

1) la durata media del processo di tunnelling TT U n(0, d) non dipende dalla profondita d 
della barriera ("effetto Hartman"); 

2) la quantita rT un (0, d) decresce all'aumentare dell'energia, come nel caso dei tempo di 
fase estrapolato; 

3) II valore della durata del tempo di penetrazione aumenta rapidamente soprattutto nella 
parte iniziale della barriera, in prossimita cioe del punto x = 0; 
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4) Tp en (0, x) tende a un valore di saturazione della parte finale della barriera, cioe per 
x — ► d. 

Per quanto riguarda gli effetti riportati nei punti 3) e 4), secondo gli autori, sarebbero 
causati dall'interferenza tra le onde iniziali che penetrano nella barriera e le onde che tornano 
indietro (sempre dentro la barriera) e la cui sovrapposizione produce i flussi J+ e J_ . Si vedano 
anche le illustrazioni (in particolare Fig.3, pag.351) in ref.[6]. 

Infatti, essendo nella parte iniziale della barriera il pacchetto d'onda di ritorno abbastanza 
grande, cio che accade e che il suo flusso estingue sostanzialmente buona perte di quelle- entrante 
del pacchetto incidente. 

AH'aumentare di x, pero, il pacchetto d'onde di ritorno va a zero piu velocemente di 
quello entrante: cio provocherebbe un aumento della durata media del processo di penetrazione 
Tp en (0, x), facendo crescere velocemente tale durata, soprattutto nella zona iniziale della barriera. 
Nella regione finale della barriera il suo aumento svanisce rapidamente e si ha l'effeto inverse 

Infine, in connessione con i grafici di tr^(x, x) in funzione di x, mostrati in Fig.I-9, 
osserviamo che: 

5) la durata media di riflessione tr(0, 0) = trS~(0, 0) non dipende dall'ampiezza d della 
barriera; 

6) in corrispondenza della regione della barriera tra e, circa, 0.6 d il valore di tr^~(0, x) 
e quasi costante; 

7) il suo valore aumenta con x solamente nella regione della barriera vicino ad x = d, 
anche se come sottolineano gli autori i calcoli relativi a tale regione non sono molto accurati a 
causa del fatto che la quantita Jf^ J-(x,t)dt assume valori molto piccoli). 

Notiamo quindi che il punto 5), previsto, come il punto 1) per particelle quasi- 
monocromatiche da Dumont e Marchioro,[32] risulta in accordo con i dati ottenuti da Steinberg 
et al.[33] per pacchetti d'onde arbitrari. Come prima anche i punti 6) e 7) possono essere spiegati 
da fenomeni di interferenza all'interno della barriera: infatti, se vicino ad x = d il pacchetto 
d'onde di ritorno iniziale e smorzato quasi totalmente dal pacchetto d'onda incidente iniziale, 
allora restera solamente una parte trascurabile della sua coda posteriore (fatta dalle componenti 
con velocita minori). 

Al diminuire di x [x — > 0) la parte non smorzata del pacchetto di ritorno sembra di- 
ventare sempre piu grande (riacquistando le componenti piu rapide) cosicche la differenza 
tr^ (0, x) — Tp^(0, x) rimane, con buona approssimazione, costante. Inoltre l'interferenza tra le 
onde incidenti e riflesse nei punti x < costituisce in effetti un fenomeno ritardante cosicche 
t-(x = 0) e piu grande di tr( x = 0)i u cne P u o spiegare i valori piu grandi di tr(x = 0, x = 0) 
rispetto a tt U ii (x = 0, x = a). 
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Parte II: ESPERIMENTI. 



II. 0) Equivalenza ottica del tunnelling. 

Come gia accennato, varie conferme sperimentali dell' effetto Hartman si sono recentemente avute 
in seguito ad una serie di misure effettuate a Colonia[34], Berkeley [35], Firenze[36], e Vienna[37]. 

Tali misure, pero, sono state eseguite sui tempi di trasmissione di microonde e fotoni, 
sfruttando la propagazione di modi evanescenti all'interno di guide d'onda sotto la frequenza di 
cut-off, nel primo caso, e la riflessione frustrata, nel secondo. 

Notiamo infatti che, malgrado in passato siano state avanzate numerose proposte di es- 
perimenti da effettuare direttamente con particelle quali ad esempio elettroni, grosse dimcolta 
erano state incontrate dal punto di vista pratico nel realizzarli, a causa soprattutto dei tempi 
molto piccoli coinvolti in tali processi di tunnelling. 

Per una giunzione Josephson, ad esempio, tali tempi risultano essere di qualche decina di 
femto-secondi, e possono scendere all'ordine del fs in altri dispositivi a stato solido. 

Nel caso di sistemi ottici, invece, tali tempi sono gia dell'ordine di qualche ps, per frequenze 
nella regione del visibile, e raggiungono il ns in alcuni degli esperimenti compiuti a Firenze ed a 
Colonia, con microonde. 

Pur senza prendere qui in rassegna tali esperimenti, oramai noti e famosi[38], ed i risul- 
tati da questi ottenuti, occupiamoci dell 'equivalenza tra trasmissione di modi elettromagnetici 
evanescenti e tunnelling di particelle, in particolare nel caso delle guide d'onda. Consideriamo 
una particella di massa m ed energia cinetica E = h 2 k 2 /2m. Nel caso (unidimensionale) di 
attraversamento di un potenziale uniforme Vo, l'equazione di Schrddinger per tale particella 
sara: 

^ + |^- W = 0. (J/.0.0) 

Posto allora: 

K 2 = ^E-V ), (J/.0.1) 

la (II. 0.0) risulta formalmente identica all'equazione di Helmholtz per la componente 
scalare relativa al campo elettrico, o a quello magnetico, di un campo (e.m.) che si propaghi in 
un mezzo dispersivo: 

#V , . 2 



dx 

dove in questo caso: 



+ K 2 ip = 0, (II. 0.2) 



2tt 2tt 

\ m e la lunghezza d'onda all'interno del mezzo, A e la lunghezza d'onda nel vuoto, e n e l'indice 
di rifrazione del mezzo in cui il campo si propaga. 
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II confronto tra le due equazioni suggerisce la sostituzione: 



12m ,„ . 2vr 



Nel caso di una guida d'onda rettangolare di dimensioni a x b (a < b), e con pareti perfettamente 
conduttrici, sappiamo che: 



(II. 0.3) 




dove A c = 26 e la lunghezza d'onda di cut-off al di sopra della quale il termine sotto radice 
divene negativo e, di conseguenza, k immaginario. 
Poiche inoltre A = cjv = 2ttc/uj, abbiamo: 

, , / /, , \ 2 

(J/.0.4) 

con u) c = iTc/b (v c = c/2b), che ci da a sua volta la frequenza di cut-off al di sotto della quale k 
diviene immaginario. 

Notiamo che la relazione di dispersione nel caso di guida d'onda rettangolare e sorpren- 
dentemente uguale[39] (si vedano anche le "Feynman Lectures" [7]) a quella di una particella 
relativistica quando facciamo la sostituzione: 

7T L0 C mC 

b c h 
infatti, moltiplicando la (II. 0.4) per H abbiamo: 




Onde evitare confusione, precisiamo che qui il confronto lo stiamo facendo tra l'equazione 
di Helmholtz, che e relativistica e classica (perche ricavata dalle equazioni di Maxwell), e 
l'equazione di Schrodinger, quantistica e non relativistica. 

Proprio questo e, inoltre, uno dei vantagi di usare campi elettromagnetici anziche parti- 
celle. Torniamo, dunque, al suddetto confronto tra la (II. 0.0) e la (II. 0.2). Differenziando la 
(II.0.1) e la (II.0.4), otteniamo: 

particella _ _ ^ 

gruppo dK m 

c.em = ^ = C 1 K 
grUPP ° dK UJ 

Attraverso la sostituzione: 

h c 2 c . . 
►- = «—» ^-0-6 

m LO ITTU 

possiamo sempre adattare i risultati elettromagnetici, nel caso di trasmissione di microonde 
in una guida d'onda, a quello quantistico, del moto unidimensionale di una particella in un 
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potenziale uniformed 9 . In entrambi i casi infatti le soluzioni delle due equazioni saranno date 
da combinazioni lineari delle funzioni d'onda: 

${x,t)=e* KX e iut 

Inoltre, quando l'energia della particella risulta minore di Vo, o quando nel caso elettro- 
magnetico la pulsazione u diviene minore di u c , k diviene immaginario e le funzioni d'onda 
divengono degli esponenziali decrescenti* 10 della forma e~\ K \ x (onde evanescenti) . 

Notiamo dunque che il campo (o equivalentemente la funzione d'onda della particella) 
penetra ugualmente all'interno della zona sotto cut-off(o classicamente proibita, per le particelle) 
almeno per una distanza dell'ordine di | /«| 1 . 

Ovviamente resta il fatto che, nonostante le analogie formali, vi sono delle differenze fisiche 
tra il tunnelling di elettroni e la propagazione guidata di microonde sotto cut-off. Infatti, come 
gia detto, l'equazione di Helmholtz (per le onde) e l'equazione di Schrddinger (per elettroni) 
sono si, sostanzialmente, la stessa equazione, ma mentre nel primo caso cio che si propaga e 
effettivamente una componente del campo (e quindi e il campo stesso), nel secondo caso si 
tratta della funzione d'onda della particella. 

Tuttavia, trattandosi in ambo i casi dell'evoluzione di pacchetti d'onda, nulla ci vieta di 
interpretare i risultati degli esperimenti su microonde come delle vere e proprie "simulazioni 
fisiche" del tunnelling di elettroni, anche perche il fatto che i risultati di tali simulazioni [34 — 38] 
riproducano abbastanza bene le previsioni quantistiche conferma Pequivalenza tra i due casi. 

Phi sottile e la circostanza che nel caso dipendente dal tempo le equazioni di Schroedinger 
(in presenza di barriera) e l'equazione di Helmholtz (per onde elettromagnetiche in guida d'onda) 
non sono piu matematicamente equivalenti, dato che la derivata temporale e del primo ordine 
in un caso, e del secondo ordine nell'altro caso. Ciononostante, si pud far vedere[40] che esse 
ammettono ancora classi di soluzioni analoghe, differenti solo per le loro proprieta di spreading. 

II. 1) Misure del gruppo di Colonia. 

Riportiamo ora i risultati di alcune delle misure eseguite a Colonia dal gruppo di Nimtz, 
Enders, et al.[34], effettuate sulla trasmissione di microonde all'interno di una guida d'onda, al 
di sotto della frequenza di cut-off. 

Figura II- 1. 

Didascalia della Fig.II-1.: 
{Configurazioni sperimentali del gruppo di Colonia.} 



°# 9 Notiamo che la (II. 0.6) cquivalc alia relazionc fiui — mc 2 

°# 10 Lc soluzioni del tipo e' K ' a: , per motivi fisici, si riferiscono al set-up sperimentale opposto. 
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Phi precisamente, sono state utilizzate due guide d'onda rettangolari (ved. Fig.II.l), 
la prima delle quali, di dimensioni 10.16x22.86 mm 2 (banda X), e stata sempre utilizzata a 
frequenze sempre superiori a quelle di cut-off, di 6.56 GHz. 

A meta di questa era inserita una seconda guida d'onda, di dimensioni 7.90x15.80 mm 2 
(banda Ku), la cui frequenza di cut-off era di 9.49 GHz. 

Onde evitare eventuali errori sistematici provenienti dal restringimento dovuto alia con- 
nessione delle due guide, o all'interferenza tra la linea di "iniezione" e di ricezione del segnale con 
la guida d'onda di banda X, gli autori hanno utilizzato una speciale tecnica di calibrazione (cali- 
brazione LRM). [Piu recentemente e apparso un nuovo lavoro degli stessi autori in cui l'apparato 
sperimentale e stato leggermente variato; infatti, per evitare eventuali problemi dovuti alia stroz- 
zatura nel passaggio tra le due guide di banda X e Ku, anziche usare una seconda guida per 
creare la zona di evanescenza, sono state inserite all'interno di questa degli strati di materiale 
con indice di rifrazione diverso, sfuttando cosf la rinessione frustrata, come nel caso degli es- 
perimenti con fotoni: i risultati di tali nuove misure concordano perfettamente con quelli qui 
riportati.] 

Figura II-2. 

Didascalia della Fig.II-2: 
{Ritardi di fase sperimentali del coefficiente di trasmissione totale T(y), in funzione della fre- 
quenza, per quattro diverse lunghezze della guida d'onda centrale di banda Ku.} 

Una prima serie di misure e stata fatta nel dominio di frequenza, misurando cioe il 
coefficiente di trasmissione totale T(y), al variare della frequenza portante di un certo inpulso 
(che considereremo piu avanti). 

La quantita T(u) include quindi, in questo caso, anche la variazione di fase causata 
dall'attraversamento della guida di banda Ku. 

In effetti, le misure della fase e dell'ampiezza sono state eseguite separatamente, anche se 
contemporaneamente. 

Le misure sono state effettuate nelPintervailo di frequenze: 8.2 — 9.2 GHz, appena 300 
MHz al di sotto della frequenza di cut-off. Tale intervallo e stato diviso in 801 punti di misura: 
la differenza di frequenza tra due punti di misura successivi essendo di 1.25 MHz. 

In Fig. II. 2 sono riportati i ritardi di fase ottenuti per quattro differenti lunghezze della 
guida d'onda centrale; piu esattamente L = 40, 60, 80, 100 mm. 

E piu che evidente che i suddetti ritardi di fase sono uguali per tutte e quattro le 
lunghezze con un'accuratezza di ±1 grado. Solo alle frequenze piu piccole appaiono delle devi- 
azioni statistiche maggiori. 

Figura II-3. 

Didascalia della Figura II-3: 
{In figura sono riportati l'impulso di riferimento (L=0 — linea punteggiata) ; l'impulso trasmesso 
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attraverso una "barriera" di 100 mm (linea continua), con un ritardo di circa 130 ps; Pimpulso 
trasmesso lungo la stessa distanza di 100 mm, ma attraverso una guda d'onda di banda X (linea 
tratteggiata).} 

Va pero notato che a basse frequenze l'attenuazione del segnale e tale da far scen- 
dere la precisione degli strumenti di misura. Inoltre, anche nel vuoto un'onda elettromagnetica 
del genere di quelle usate presenta normalmente una variazione di fase dell'ordine di 1 grado, 
dopo aver viaggiato per una distanza di appena 0.1 mm. L'indipendenza dello sfasamento dalla 
lunghezza delle guide intermedie dimostra che questo e causato essenzialmente solo dalle con- 
dizioni al contorno sulle due facce (d'ingresso e di uscita) della guida d'onda di banda Ku. 

Figura II-4. 

Didascalia della Fig.II-4: 
{Impulso di riferimento (L=0), linea continua; impulsi di tutte e quattro le misure, che sono 
circa uguali, linee punteggiate; impulso che attraversa sempre i 100 mm ma in una guida di 
banda x al di sopra quindi della frequenza di cut-off, linea tratteggiata.} 

Le stesse misure sono state riportate, dagli autori, nel dominio temporale attraverso 
un'integrazione di Fourier. Supponiamo infatti di avere un segnale F(t), dato da : 

/'OO 
-oo 

dove f(u) e la trasformata inversa di Fourier del segnale F(t). Dopo che il segnale avra viaggiato 
attraverso una certa regione di spazio, la sua nuova forma, nel dominio del tempo, sara: 

/oo 
f(v)T(v)e i2nut dt. (J/.1.0) 
-oo 

Dunque, una volta misurata T(u), e conosciuta la forma del segnale iniziale, e possibile 
ricavare direttamente le misure dei tempi di arrivo del segnale. Gli autori adottarono come 
forma del segnale quella di una distribuzione di Gauss, supponendo, inoltre, che fosse possibile 
limitare Pintegrale (II. 1.0) all'mtervallo di frequenze (^1,^2), dove v\ e v<i sono le frequenze di 
cut-offdelle due guide d'onda. 

In Figg.II.3 e II. 4 sono riportati, rinormalizzati, i risultati ottenuti trasformando le 
precedenti misure nel dominio del tempo (ovvero, in pratica, le F'(t)). Notiamo come la de- 
formazione degli impulsi attraverso la regione della barriera sia trascurabile: non c'e quindi 
assolutamente bisogno di distinguere tra picco e centra di massa dell'impulso. 

Figura II-5. 
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Didascalia della Fig.II-5: 
{Inviluppi delle potenze dei segnali trasmessi per due frequenze portanti: a) 8.44 GHz, b) 8.658 
GHz. Le figure continue rappresentano gli inviluppi dei segnali prima della barriera, attenuati 
di 40 dB.} 

Gli stessi autori eseguono poi una seconda serie di misure direttamente sui tempi di 
attraversamento mediante un analizzatore di tempi. In questo caso, il segnale viene modulato 
in ampiezza creando un treno di impulsi, con frequenza portante v$. La durata di ogni impulso 
e di un /xs, mentre la frequenza di ripetizione di 10 kHz. 

L'analisi di Fourier del segnale comprende, in questo caso, uno spettro infinito di frequenze, 
distribuite, comunque, intorno alia frequenza portante uq. Tale frequenza va dunque scelta in 
modo da evitare che Pimpulso trasmesso sia prevalentemente composto da frequenze al di sopra 
di quella di cut-off. 

Siccome pero queste frequenze non vengono attenuate, generano distorsioni abbastanza 
vistose del segnale. Secondo gli autori, e abbastanza semplice, verificare quando cio accade, e 
mettersi in condizioni da poterlo evitare (cf. la prima delle refs.[34]). 

Le misure sono state fatte in questo caso sui tempi di arrivo del fronte d'onda di ogni 
singolo impulso. 

In Figg.II.5a e II. 5b, sono rappresentati (sia in scala naturale, a sinistra, sia in scala 
logaritmica, a destra), in funzione di t, gli inviluppi delle intensita di due impulsi, con frequenze 
portanti rispettivamente di 8.644 e 8.658 GHz, per uno spessore della guida centrale di 60 mm. 
Le linee continue invece rappresentano lo stesso inviluppo per il segnale di riferimento (L = 0). 
Quest 'ultimo veniva fatto passare attraverso un attenuatore di potenza a scalino, che ne riduceva 
l'ampizza di 40 dB. L'errore temporale introdotto dall'attenuatore e comunque inferiore a 25 ps, 
indipendentemente dalla frequenza. 

Come si puo vedere, anche in questo caso, il tempo di attraversamento risulta indipendente 
dallo spessore della guida "subcritica" centrale. 

Misure analogue sono state fatte, sempre dagli stessi autori, anche per altri valori della 
lunghezza della guida centrale, ottenendo risultati simili a quelli delle due figure qui riportate. 

Figura II-6. 

Didascalia della Fig.II-6: 
{Rappresentazione schematica della configurazione sperimentale nel caso della doppia barriera.} 

II caso piu interessante studiato dallo stesso gruppo, infine, e stato quello 
dell'attraversamento di una doppia barriera (cioe di due tratti di guida d'onda subcritica). 
Riportiamo qui i risultati sperimentali ottenuti in questo caso, in quanto i tempi di attraversa- 
mento misurati, per bande di frequenze lontane dalle frequenze di risonanza, non solo risultano 
essere indipendenti dallo spessore delle due barriere, ma sono indipendenti anche dalla distanza 
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tra queste due. 



Figura II-7. 
Didascalia della Fig.II-7: 

{La figura a destra riporta le intensita trasmesse di un pacchetto d'onda gaussiano in funzione 
del tempo nel caso di doppia barriera, per le tre configurazioni [a), b), c)] adottate. Notare 
l'equivalenza delle configurazioni a) e b). Sulla sinistra invece e riportato il coefficiente di 
trasmissione delle configurazioni a) e b). Per la configurazione b) si vedono apparire due fre- 
quenze di risonanza.} 

In Fig.II.6, possiamo vedere le configurazioni sperimentali adottate: 

a) dapprima, le due guide di banda Ku, entrambe di lunghezza 40 mm, sono affiancate 
l'una all'altra; 

b) successivamente, una terza guida di banda X e lunghezza 57 mm viene posta tra le 
due guide sottodimensionate; 

c) la guida di 57 mm viene, infine, rimossa e montata in coda alia altre, mentre le due 
guide di banda Ku vengono riaccostate. 

Con la configurazione b) appaiono due frequenze di risonanza, rispettivamente a 6.9 e 
7.6 GHz. In prossimita di tali risonanze, i tempi di tunnelling tunnelling misurati dagli autori 
tendono al valore inverso della larghezza della risonanza. 

Piu importante e la misura del tempo di attraversamento di un pacchetto gaussiano, 
stretto in k (Ak = 0.01 A;), alia frequenza non risonante di 7.3 GHz. I risultati di tale misura 
sono riportati in Fig.II.7 (grafico a sinistra). Come si puo vedere, il tempo di attraversamento 
del segnale e uguale sia nel caso a), senza "buca di potenziale", sia nel caso b) con una guida di 
banda X di 57 mm posta tra le due guide di banda Ku. 

Nel caso c) il tempo di attraversamento dei 57 mm risulta di 400 ps in accordo con la 
velocita di gruppo alPinterno della guida. 

Dalla seconda delle due Figure II. 7, notiamo che il coefficiente di trasmissione, lontano 
dalle frequenze di risonanza, assume all'incirca gli stessi valori, sia nel caso a), sia nel caso b). 

I due autori passano quindi a considerare — teoricamente — il caso della misura con 
la configurazione b), ma con una guida intermedia lunga 1 km (invece di 57 mm). In tal caso, 
naturalmente, il numero delle risonanze aumentera, diminuendo la distanza tra l'una e Paltra. 
II pacchetto d'onda usato per misurare il tunnelling non risonante va dunque stretto ad una 
larghezza di 10~ 6 volte la frequenza centrale, onde evitare effetti di interferenza. Allora, per 
una frequenza portante di 10 GHz, sarebbe possibile, ad esempio, trasmettere un segnale di 
larghezza 10 KHz. Secondo gli autori un segnale del genere basterebbe per modulare e trasmet- 
tere anche la Sinfonia n.40 di Mozart ad una distanza di 1 km ed a velocita Superluminalc. 
Gli autori si riferiscono a simulazioni attraverso computer [41], di cui danno pero pochi dettagli. 
Tali computer simulations sono state pertanto rifatte con cura presso il D.M.O. della Electric 
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Engineering Faculty della Campinas State University, Campinas, S.P., Brasile[42], con analoghi 
risultati di cui si daranno i risultati altrove. 



Figura IF8. 

Didascalia della Fig.II-8: 
{Grand dei tempi di trasmissione calcolati da Hartman, in funzione dello spessore della bar- 
riera (de), e del numero d'onda n' rinormalizzato a e (k' = k/s, essendo e il numero d'onda 
corrispondente all'altezza della barriera). I punti segnati rappresentano approssimativamente i 
dati sperimentali per un pacchetto elettromagnetico gaussiano, con frequenza centrale 8.7 GHz, 
nel caso di spessori della guida centrale (frequenza di cut-off9.49 GHz) di 10, 40, 60, 80, 100 
mm.} 

Riportiamo, infine, i grafici dei tempi di tunnelling di una barriera di potenziale di 
spessore d ed altezza H 2 s 2 /2m, calcolati da Hartman. Nello stesso grafico sono riportati i rel- 
ativi tempi equivalenti. I punti segnati rappresentano approssimativamente i dati sperimentali 
per un pacchetto elettromagnetico gaussiano, con frequenza centrale 8.7 GHz, nel caso di spes- 
sori della guida centrale (frequenza di cut- off) .49 GHz) di 10, 40, 60, 80, 100 mm. Notiamo 
l'impressionante precisione dei dati sperimentali. 

II. 2) Misure del gruppo di Berkeley. 

Prendiamo ora in considerazione i famosi esperimenti condotti a Berkeley dal gruppo di 
Chiao, Steinberg e Kwiat. 

In questo caso le misure sono state svolte sfruttando la band-gap dovuta alle riflessioni 
frustrate su un materiale multistrato. Anche in questo caso, comunque, vale l'analogia tra 
l'equazione di Schrodinger e quella di Helmholtz, e possiamo estendere i risultati delle misure 
al caso di barriere di potenziale quantistiche. Quella che cambia sara, al solito, la relazione di 
dispersione del mezzo e, come nel caso delle guide d'onda, anche qui, in un certo intervallo di 
frequenze, il vettore d'onda k diviene immaginario, e appaiono i modi evanescenti. 

L'uso di fotoni, pero, presenta alcuni vantaggi, rispetto a quello delle microonde: per 
esempio, il fotone dimostra un comportamento "individuale" diverso da quello presentato da 
un segnale elettromagnetico qualsiasi nel campo delle microonde. Inoltre anche la dispersione e 
molto piii ridotta, quasi come nel caso di una particella. 

Figura II-9. 

Didascalia della Fig.II-9: 
{Configurazione sperimentale di Berkeley.} 
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In Figura II. 9 e schematizzato l'apparato sperimentale usato in questo caso. 
Una coppia di elettroni viene generata attraverso un processo di riconversione spontanea 
provocato da un laser all'interno di un cristallo (KDP) di potassio-di-idrogeno-fosfato. I due 
elettroni sono generati contemporaneamente e mentre uno dei due viaggia liberamente in aria, 
l'altro viene fatto incidere su un supporto di materiale rifrangente, meta del quale e ricoperto 
da un multistrato composto di due sostanze con indice di rifrazione diverse Ogni strato ha uno 
spessore pari a un quarto della lunghezza d'onda dei fotoni. Lo spessore totale e di 1.1/xm, e 
corrisponde ad un tempo di attraversamento d/c = 3.6 fs. La band-gap dovuta al multistrato si 
estende per un'intervallo di lunghezze d'onda che va da 600 a 800 nm, con un minimo dell'1% 
a 692 nm nella trasmissione. La faccia opposta e invece ricoperta di materiale antiriflessione. 
Spostando il supporto e possibile introdurre, o no, una zona di evanescenza sul percorso del 
fotone. I due fotoni vengono poi fatti ri-incrociare, attraverso alcuni specchi, in modo da farli 
passare contemporaneamente per un divisore di fascio al 50 %. Sulle uscite del divisore di fascio 
sono posti due rivelatori, che registrano una coincidenza se i due fotoni arrivano entro 500 ps 
l'uno dall'altro. 

Quando i due elettroni arrivano simultaneamente all'interno del divisore di fascio, 
l'interferenza tra i due ha un effetto distruttivo, annullando cosf la coincidenza. Aumentando o 
diminuendo, attraverso un prisma mobile, il cammino di uno dei due elettroni e possibile sta- 
bilire quando i due attraversano contemporaneamente il divisore di fascio, o quale dei due arriva 
prima. 

Introdotto quindi lo spessore di multistrato sul cammino di uno dei due elettroni, e pos- 
sibile misurarne, in questo modo, l'anticipo o il ritardo che questo ha rispetto all'altro. 

Figura 11-10. 

Didascalia della Fig. II- 10: 
{Profili dei conteggi di coincidenze con e senza la barriera riflettente (rispettivamente curva piu 
in alto e curva piu in basso). Come si puo vedere la curva phi alta mostra un anticipo di circa 
1.1 (± 0.3) fs rispetto all'altra.} 

In Figura 11.10 e riportato il risultato ottenuto dai tre autori sia nel caso di presenza 
della "barriera" , sia in sua assenza. 

Si puo dunque vedere che nel caso di presenza del multistrato, i fotoni che lo attraversano 
arrivano, in media, circa 1.1 fs prima degli altri. L'errore considerato dagli autori e di 0.3 fs. 
Cio dimostra che che il tempo impiegato da ogni singolo fotone ad attraversare la barriera e 
di circa 7 deviazioni standard piu piccolo del tempo che lo stesso elettrone impiegherebbe per 
attravesare lo spazio da questa occupato, ma in sua assenza, alia velocita c. 

Gli autori fanno notare inoltre che, in questo caso, poiche per quasi tutto l'intervallo 
della band-gap il coefficiente di trasmissione varia molto lentamente, mantenendosi quasi 
costante al variare dell'energia, non e possibile giustificare tale evento come apparente e dovuto 
all'attraversamento, preferenzialmente, di fotoni piu energetici. Inoltre notiamo che tutti i fotoni 
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nell'aria viaggiano circa alia stessa velocita, molto prossima a quella della luce, indipendente- 
mente dalla loro energia. 

Gli stessi autori, pero, sostengono la tesi opposta, in un articolo divulgativo apparso su 
Scentific American nell'Agosto del 1993. 

II. 3) Misure del gruppo di Firenze. 

Pur operando anch'esso con microonde, il gruppo di Firenze riprende sia dal punto di 
vista delPanalisi teorica, che nella procedura sperimentale un'idea di Brillouin[*41], sviluppata 
in seguito da Stevens [*42, *43]. Secondo questi infatti, in un mezzo dispersivo, va considerata 
come velocita di propagazione di un segnale elettromagnetico la velocita di propagazione del 
fronte d'onda molto netto di un segnale a scalino, e non la sua velocita di gruppo. 

Notiamo che, anche se un segnale del genere puo essere anch'esso soluzione dell'equazione 
di Helmholtz, non esiste pero un equivalente quantistico. La propagazione di un segnale del 
genere e, poi, molto piu dispersiva di quella di altri tipi di segnali, ed e accompagnata sempre 
dall'apparire, dopo un certo intervallo di tempo, di alcuni precursori, formati dalle frequenze 
piu alte. 

Figura 11-11. 

Didascalia della Fig.II-11: 
{Ritardi temporali in funzione della frequenza per una barriera di 15 cm di lunghezza. In figura 
sono pure riportate le curve teoriche di r s , e r^.} 

Tralasciamo dunque le analisi teoriche suH'argomento, ma prendiamo ugualmente in 
esame i loro dati sperimentali, anche perche in essi ci troviamo un confronto diretto con alcuni 
dei tempi precedentemente definiti nel primo capitolo. 

La procedura usata assomigia abbastanza a quella del gruppo di Enders e Nimtz nella 
prima delle loro misure riportate. Le guide d'onda sono in questo caso rispettivamente: di 
banda X quella esterna e di banda P quella interna (per le dimensioni vedere ref. [*44]). In 
una prima serie di misure, un segnale a scalino, nell'intervallo di frequenze 9.4, 9.7 GHz viene 
fornito da un klinstron, mentre il segnale trasmesso viene mandato ad un oscilloscopio ad alta 
risoluzione per la misura dei ritardi di fase e dell'attenuazione del segnale. 

I risultati dei ritardi di fase sono corretti dagli autori mediante la sottrazione del ritardo 
corrispondente alia misura fatta senza la guida di banda P. 

In Figura 11.11 sono riportati i valori dei ritardi di fase misurati in funzione delle fre- 
quenze, insieme con le curve teoriche di t s , t% e . Vediamo, che anche se nessuna delle curve 
riproduce esattamente i dati sperimentali, quella del phase time sembra si avvicini abbastanza 
ai dati. 
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Figura 11-12. 



Didascalia della Fig.II-12: 
{Come per la Figura 11.11, ma per L = 10 cm.} 

La Figura 11.12 si riferisce, invece, al caso di L = 10 cm. In questo caso pero i punti 
sperimentali ricadono tra la curva relativa ar| e quella del phase time, mentre r s rimane in- 
adeguato. Gli autori attribuiscono tale comportamento alia giunzione tra le due guide, infatti, a 
quanto pare, non viene usata nessuna tecnica di calibrazione del sistema. L'effetto aumenterebbe, 
secondo gli autori, al diminuire della lunghezza della barriera. 

II. 4) Conclusioni. 

Come abbiamo visto, il problema posto all'inizio, di quanto tempo impieghi una particella ad 
attraversare una zona classicamente proibita, benche sostanzialmente risolto, rimane ancora 
controverso. 

Se da un lato infatti esiste una evidenza sperimentale di gran lunga a favore della 
definizione di phase time, dall'altro questa stessa defmizione, e gli stessi risultati sperimen- 
tali, portano a dover accettare l'insorgere di velocita di gruppo maggiori di quella della luce nel 
vuoto, a causa di quello che Olkhovsky e Recami[5] hanno chiamato "effetto Hartman". 

Va comunque notato che, l'apparizione di velocita di gruppo Superluminali o addirittura 
negative, nell'ottica classica, non e un fenomeno nuovo, e fu anche questo studiato ed osservato 
sperimentalmente in lavori[45*, 47*] che solo di recente stanno richiamando l'attenzione che 
meritavano. Infatti, come sappiamo, in un mezzo dispersivo lineare, non assorbente: 

du c 
s dk n(uj) + un'{uj) ' 

e, in regioni forte dispersione anomala (in prossimita di risonanze), la velocita di gruppo puo 
eccedere c o, come gia detto, divenire negativa. E comunque dimostrato, pero, che in questi casi 
la velocita di gruppo perde di validita fisica e nessun segnale puo in effetti essere trasportato 
dal mezzo a velocita maggiore di c [si vedano anche il Cap. 7 del Classical Electrodynamics del 
Jackson, o il Cap. 3 del Sommerfeld]. 

In tal caso, inoltre, il fenomeno e facilmente spiegato attraverso un reshaping dell'impulso, 
che subisce un'attenuazione delle componenti meno energetiche e piu lente. II segnale uscente 
sembra quindi viaggiare a velocita apparentemente maggiori di quella a cui hanno realmente 
viaggiato le sue componenti phi energetiche. 

Nel tunnelling, pero, abbiamo visto che (sia per le particelle, sia nel caso elettromagnetico) 
possiamo sempre metterci in condizioni da evitare effetti dovuti ad eventuali reshaping, o alia 
trasmisiione di particelle gia inizialmente piii veloci. 
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Cio e anche confermato sperimentalmente dal fatto che, come abbiamo visto soprattutto 
nelle misure di Enders e Nimtz, la forma dei segnali rimane inalterata. 

Notiamo pero che il processo di attraversamento della barriera, da parte di una particella, e 
comunque un processo statistico, nel senso che non possiamo mai sapere a priori quale particella 
attraversa la barriera. Inoltre, l'apparire di velocita Superluminali si ha proprio quando la 
probability associata a tale evento e molto bassa. 

Secondo alcuni autori, quindi, una giustificazione della violazione del principio di causalita 
potrebbe essere che tale fenomeno, in quanto non controllabile, non puo essere usato per trasmet- 
tere alcuna informazione. 

In quest'ottica si inserisce il tipo di giustificazione che portano Steinberg et al.[10]. Sec- 
ondo questi ultimi, nei processi di tunnelling in questione (quelli Superluminali), i relativi tempi 
di attraversamento andrebbero visti come "valori deboli" della misura. II concetto di valore 
debole, o misura debole, e stato introdotto nel 1988 da Aharonov, Albert e Vaidman[49*, 50*], 
partendo dalla teoria "classica" della misura di von Neumann [51*]. 

Secondo i due autori quando facciamo una misura, su di un sottoinsieme con bassa prob- 
abilita associata (come quelo delle particelle trasmesse), e questo sottoinsieme proviene da un 
insieme su cui e stata effettuata una "misura debole" (una misura cioe con una grossa indermi- 
nazione, che lasci lo stato del sistema quasi imperturbato), e possibile ottenere come risultato 
della misura sul sottoinsieme, un valore completamente diverso da tutti gli autovalori accessibili 
al sistema. Tale valore non sarebbe comunque un valore realmente assunto dal sistema. 

Secondo altri autori, invece, le velocita Super-luminali legate al tunnelling sarebbero pro- 
prio reali e quello che va reinterpretato e il principio di causalita. [46*] A quest'ultimo proposito, 
ricordiamo che la Relativita Speciale puo essere estesa — senza variarne i consueti Postulati — 
in modo da inglobare i moti Superluminali; la "Relativita Estesa" [47*] , in altre parole, include 
i tachioni senza sostanzialmente violare la Relativita di Einstein, ma solo estendendola al nuovo 
dominio di velocita. In particolare, si possono risolvere le questioni causali.[48*] 

E anche utile ricordare che la stessa Relativita Estesa prevede per oggetti Superluminali, 
sulla base di semplici cosiderazioni geometrico-classiche, la transizione da velocita di gruppo 
positive a velocita negative (quando si "oltrepassi" la situazione di velocita infinita): si vedano 
le refs.[47*]. Dato che a velocita negative corrispondono tempi di transito negativi[29*], cio 
diviene interessante alia luce di esperimenti meno recenti[42*] e phi recenti [34-38]. 
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